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Abstract 

We study approximate solutions to the time-dependent Schrodinger 
equation iedt'>ptix)/dt — H{x, —ie\7 x) "fptix) with the Hamiltonian given 
as the Weyl quantization of the symbol H{q,p) taking values in the 
space of bounded operators on the Hilbert space Hf of fast "inter- 
nal" degrees of freedom. By assumption H{q,p) has an isolated energy 
band. Using a method of Nenciu and Sordoni [NeSo] we prove that in- 
terband transitions are suppressed to any order in e. As a consequence, 
associated to that energy band there exists a subspace of L^(R'^,7if) 
almost invariant under the unitary time evolution. We develop a sys- 
tematic perturbation scheme for the computation of effective Hamil- 
tonians which govern approximately the intraband time evolution. As 
examples for the general perturbation scheme we discuss the Dirac and 
Born-Oppenheimer type Hamiltonians and we reconsider also the time- 
adiabatic theory. 
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1 Introduction 

Quantum theory has the remarkable feature that certain dynamical degrees 
of freedom may become "slaved" and thus lose their autonomous status. 
The origin of this phenomenon is a separation, both in space and time, into 
slow and fast degrees of freedom. The fast modes quickly adapt to the slow 
modes which in turn are governed by a suitable effective Hamiltonian. This 
mechanism is called adiabatic decoupling. 

As paradigm we mention the motion of nuclei. The electronic degrees 
of freedom rapidly adjust to the state of lowest energy at given positions 
of the nuclei and the electronic energy band serves as effective potential in 
the Hamiltonian for the nuclei. This Born-Oppenheimer approximation is 
the basis for the dynamics of molecules and, as a consequence, also for the 
microscopic theory of classical fluids. There are many other examples of a 
similar structure. A very widely studied case are electrons moving in the 
periodic crystal potential, which defines the short scale. The envelope of the 
electronic wave function is governed by an effective Hamiltonian obtained 
from the Peierls substitution, in which the band energy is taken as effective 
kinetic energy. For an electron coupled to the quantized radiation field the 
photons are the fast degrees of freedom and the dynamics of the electron is 
governed by an effective Hamiltonian accounting for spin precession. These 
and other systems have been studied extensively by model specific approx- 
imate methods without realizing that they share a common structure. In 
fact, all the examples given can be molded into the generic form 

d 

i£—ipt{x) = H{x, -ieVx)ipt{x) ■ (1) 

Here H{q,p) is an operator- valued function on the classical phase space 
r = M^'^ of the slow degrees of freedom with q a position like and p a 
momentum like variable. H(q,p) is self-adjoint and acts on the Hilbert space 
Hi of "internal" fast degrees of freedom. After quantization H{x, —ieVx) 
becomes the Hamiltonian of our system. To properly define it one has to 
specify an ordering of the operators x, —ieVx, for which we will adopt the 
Weyl quantization rule as to be explained in full detail in the Appendix, ipt 
is a wave function on R'^ with values in H{. Thus the quantum mechanical 
Hilbert space of states is H = L'^{R'^,ni) = L'^{R'^) ® H{. Finally e is a 
dimensionless parameter which controls the scale separation, e -C 1. 

Examples will be given in due course and we only remark that, in general, 
(1) is already the result of a proper identification of the slow degrees of 
freedom. For example, in nonrelativistic QED q stands for the position of 
the electron, whereas p is the total momentum of the electron and photons. 
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For an electron in a periodic potential (1) is the Hamiltonian in the crystal 
momentum basis with x standing for the Bloch momentum k in the first 
Brillouin zone, see [PST2]. 

Our goal is to construct approximate solutions to (1), a task which we 
divide into four steps. 

(i) Almost invariant subspace. The adiabatic decoupling can be traced 
to a spectral property of H{q,p). It is assumed, and it can be proved in 
many particular cases, that for each q,p the spectrum, a{q,p), of H{q,p) 
can be decomposed into a "relevant" part cr^{q,p) and a remainder in such 
a way that the relevant energy band {{q,p,X) S M^'^"'"^ : A € cri.(q,p)} is 
separated from the remainder {{q,p,X) € R^'^^^ : A G <Tr{q,py} by a gap, 
i.e. by a corridor of finite width. In many cases of interest the relevant 
part of the spectrum consists of a single, possibly degenerate eigenvalue and 
the relevant energy band is the graph of a smooth function. Of course, a 
given H(q,p) could have several such relevant energy bands and we suppose 
one of them to be singled out, e.g., through the initial condition. To the 
relevant energy band of H one associates a subspace UTi of Ti with the 
property that if tpQ G 117^, then tpt G IITY up to an error which is smaller 
than any power of e. For this reason Il7i is called an almost invariant 
subspace. If H denotes the Weyl quantization of H, then [-ff, H] = 0{e°°), 
where 11 is the orthogonal projector onto the closed subspace HH. Of course, 
in practice, approximations to 11 are constructed, for which at order n the 
error is 0{e"^~^^). 

(ii) Reference Hilbert space. Clearly, on IlTi the dynamics is generated 
by the diagonal Hamiltonian Ilifll, up to 0(e°°). While properly defined, 
for a further analysis of the effective dynamics inside UH a representation in 
terms of adapted coordinates for the slow degrees of freedom is more pow- 
erful. We call the corresponding Hilbert space /C = L2(M°', /Cf) the reference 
space. By definition /C is independent of e. The next task is thus to construct 
a linear map U : IIH ]C, which is unitary in the sense that U*U = 11 and 
UU* = Ik. 

(iii) Effective Hamiltonian. The effective Hamiltonian is defined through 
h = UHHUU* as operator on /C. It is unitarily equivalent to the diagonal 
Hamiltonian, but acts on a simpler space. While h is still a rather abstract 
object, it can be expanded in powers of e, which is the adiabatic perturbation 
theory of the title. Already the lowest order approximations provide a wealth 
of information on the motion of the slow degrees of freedom. 

(iv) Semiclassical limit. In many applications of interest the effective 
Hamiltonian is of a form which allows for a semiclassical analysis. The 
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simplest situation is a relevant band which consist only of a single, possibly 
degenerate, eigenvalue. Then the effective Hamiltonian has a scalar principle 
symbol and the semiclassical analysis is straightforward. 

We are certainly not the first ones to investigate approximate solutions to 
(1) and we have to explain which parts of the program outlined above have 
been achieved before and which parts are our novel contribution. In addition 
there have been other approaches to (1) on which we briefly comment below. 

To our knowledge, the notion "almost invariant subspace" was first coined 

by Ncnciu [NC2] in the context of gauge invariant perturbation theory. In 
the context of spacc-adiabatic problems Brummelhuis and Nourrigat [BrNo] 
construct 11 for the particular case of the Dirac equation and Martinez and 
Sordoni [MaSo] based on [So] consider Born-Oppenheimer type Hamiltoni- 
ans. The general scheme for the construction of 11 is sketched in Nenciu and 
Sordoni [NeSo] and applied to the matrix-valued Klein-Gordon equation. 
Our construction is based on the one in [NeSo] , but differs in a few technical 
details. 

As in the case of 11 we construct the unitary U in several steps. First 
we compute order by order the formal symbol of U, which, after quantiza- 
tion, gives rise to an "almost unitary" operator. Finally the almost unitary is 
modified to yield a true unitary exactly intertwining IlTi. and /C. Our method 
is specifically designed to deal also with problems as the Dirac equation and 
the Bloch electron with external magnetic fields, where the projector 11 has 
no limit for £ ^ 0, see Remark 3.2. While the specific application and the 
proof are new, the general idea to construct a pscudodifferential operator 
which is almost unitary and diagonalizes a given pseudodifferential opera- 
tor has a long tradition, [Ni] Section 7 and references therein, [Ta],[IIeSj]. 
The method of successive diagonalization is also prominent in the physics 
literature, for example [FoWo] in the derivation of the Pauli equation and 
its corrections, [Bli] for periodic Schrodinger operators, [BI2] for the Dirac 
equation, [LiFl, LiWe] for Born-Oppenheimer type Hamiltonians. 

Our central result is the expansion of the effective Hamiltonian h. We 
provide a scheme which is applicable in general and work out explicitly the 
expansion as h = ho + ehi + £^h2 + 0{£^). In particular, for a relevant 
band consisting of a single eigenvalue, hQ is the Peierls substitution and hi 
contains among other things the information on geometric phases. 

Since h has, in general, a matrix-valued symbol, we discuss and apply 
some results on the semiclassical limit for matrix-pseudodifferential opera- 
tors with scalar principal symbol. We include this part to make the paper 
self-contained and to demonstrate that in many cases of interest the motion 
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of the slow degrees of freedom can be approximately described through the 
appropriate classical Hamiltonian flow. 

Surely the reader will have noticed that our program makes no mention 
of the initial conditions for (1). The reason is simply that our estimates 
are uniform on UTi, respectively on /C. Physically this partial independence 
from the initial conditions is most welcome, because in general they cannot 
be controlled so easily. On the other hand, in the approach of Hagedorn and 
Joye [HaJo2] one constructs for a specific e-dependent ipQ the wave packet 
t tpf G H, which solves (1) up to or even 0(e^^/^). Hagedorn 

and Joye study Born-Oppenheimer Hamiltonians. Periodic Schrodinger op- 
erators are considered in [GRT, DGR], who employ a related multi-scale 
analysis. Whereas in the other approaches the adiabatic and the semiclas- 
sical limit are taken simultaneously, both limits are clearly separated in our 
results. We consider this an important conceptual advantage. 

To give a brief outline of our paper. In Section 2 we explain the pre- 
cise assumptions on the Hamiltonian H and construct the almost invariant 
subspace UTi. The reference space /C and the unitary intertwining it with 
IlTi. are explained Section 3. The central part of our paper is the expansion 
of the effective Hamiltonian in Section 4. In particular, we work out the 
expansion for Born-Oppenheimer Hamiltonians including order e^. We also 
show that by the Howland trick the standard time-adiabatic theory can be 
subsumed under (1). The semiclassical analysis for the effective Hamilto- 
nian is summarized in Section 5. Last but not least, in Section 6 we discuss 
the Dirac equation with slowly varying external vector potentials, since it 
is the simplest Hamiltonian for which the full generality of our approach is 
needed and yields interesting physical results. In this case TC^ = and 
the classical symbol has two two-fold degenerate energy bands, one for the 
electron and one for the positron. Thus the reference Hilbert space e.g. for 
the electron band is L^(R^,C^). The effective Hamiltonian is determined 
including order e. The principal part /iq describes the translational motion 
through the Peierls substitution and the subprincipal part hi yields the spin 
precession as governed by the BMT equation. Since the external vector 
potential appears through minimal coupling, the projection P{q,p) on the 
electron subspace depends nontrivially on both coordinates, in contrast to 
Born-Oppenheimer and time-adiabatic theory, where the relevant projection 
P{q,p) depends only on one of the canonical coordinates. We end the paper 
with some concluding remarks in Section 7 and with an Appendix reviewing 
some results on pseudodifferential calculus with operator- valued symbols. 

Acknowledgments: We are grateful to Andre Martinez and Gheorghe 
Nenciu for helpful discussions, which in part initiated our work. We thank 
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Bcrndt Thaller for his remarks on the problem of observables for the Dirac 
equation and Patrick Gerard for pointing out to us reference [Ni]. S.T. 
also thanks Raymond Brummelhuis, George Hagedorn, and Alain Joye for 
valuable comments and their interest in this work. G.P. is grateful to Fabio 
Pioli for a useful remark about fiber bundles. 



2 General setting and construction of the almost- 
invariant subspace 

Space-adiabatic perturbation theory deals with quantum systems in which it 
is possible to distinguish between fast and slow degrees of freedom. In par- 
ticular wc assume that the Hilbcrt space Ti. admits a natural decomposition 
as 7^ = L^(M'^)(8) 7^f, where L'^(M.'^) is the state space for the slow degrees of 
freedom and H{ is the state space for the fast degrees of freedom. 

As the second structural ingredient we require that the Hamiltonian is 
given as the quantization of a B(Hf) -valued function on the classical phase 
space M?*^ of the slow degrees of freedom. Hence we need to consider the 
generalization of the usual quantization rules to the case of i3(7if)-valued 
functions on M?'^. This theory is well covered in the literature, see for example 
[Ho, Fo, Iv, GMS]. Still, for the convenience of the reader, we provide a 
self-contained review of the basic results in the Appendix, where we also 
introduce the relevant notation. 

We now state the general assumptions on which the adiabatic perturba- 
tion theory will be based in the following. Let 7if be a separable Hilbert 
space, the state space for the fast degrees of freedom, and H = L^(R'^) 0H{. 
The Hamiltonian H of the full system is given as the Weyl quantization 
of a semiclassical symbol H G S'^{£). We assume that H is essentially 
self-adjoint on S. A point in the classical phase space R^*^ is denoted by 
z = {q,p) G M^"^. 

The adiabatic decoupling relies on a gap condition for the principal sym- 
bol Ho of H. 

Condition (Gap)o-. For any z G M.'^^ the spectrum a{z) of Ho{z) G B{Hi) 
contains a relevant subset aj-{z) which is uniformly separated from its com- 
plement cr{z) \ (Jj-{z) by a gap. More precisely there are two continuous 
functions -fj : M^*^ — >■ M ( j = ±) (with 7_ < 7+) such that: 

(Gl) for every z G M^*^ the spectral component ar{z) is entirely contained in 
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the interval l{z) := [7- (z), 74.(2;)]; 

(G2) the distance between a{z) \ ar{z) and the interval 1(2;) is uniformly 
bounded away from zero and increasing for large momenta, i.e. 

dist{a{z)\a,{z),Iiz))>Cg{pr; (2) 
(G3) the width of the interval 1(2;) is uniformly bounded, i.e. 

sup I7+ (2;) — 7_ (z) I < C < CXD . 

We denote the spectral projector corresponding to ar{z) by 710(2). As 
explained in the Introduction, one expects interband transitions to be sup- 
pressed for small e. To prove such a property we need either one of the 
following assumptions to be satisfied. 

Condition of increasing gap (IG)^- 

Let H be an hermitian symbol in S^{£,B{Hi)) (with p > and m > 0) 
such that the principal symbol Hq satisfies condition (Gap)(j with a = m. 

Condition of constant gap (CG). 

Let H be an hermitian symbol in (e, B{TCi)) such that the principal symbol 
Ho satisfies condition (Gap)(j with a = 0. 

Note that for the case H G Sl{e, B{Hi)) one can show that H - the Weyl 
quantization oi H - is essentially sclf-adjoint on the domain 5(]R'^,7if) C H,. 
The proof is postponed to an appendix to the space-adiabatic theorem. 

In analogy with the usual time-adiabatic theorem of quantum mechanics, 
see Section 4.4, we baptize the following result as space-adiabatic theorem. It 
establishes that there are almost invariant subspaces associated with isolated 
energy bands. In spirit the result is not new. However, to our knowledge 
it appears in this explicit form only recently in the literature. Brummel- 
huis and Nourrigat [BrNo] gave a proof for the Dirac equation, Martinez 
and Sordoni [MaSo] considered Born-Oppenheimer type Hamiltonians (cf. 
Section 4.3) based on results from [So], and Nenciu and Sordoni [NeSo] 
sketched the general scheme and applied it to a matrix-valued Klein-Gordon 
type problem. 

Theorem 2.1 (Space-adiabatic theorem). Assume either (IG)m or (CG). 
Let H he the Weyl quantization of H. Then there exists an orthogonal pro- 
jector n G B{H) such that 

[H,U]=Ooien (3) 
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and n = TT + Oo{e°°), where tt is the Weyl quantization of a semiclassical 
symbol 

TT X ^ e%j in Sp{£) 

j>0 

whose principal part tto{z) is the spectral projector of Hq{z) corresponding 
to crr(z). 

The subspace Ranll C H is thus an almost invariant subspacc for the 
dynamics generated by the Hamiltonian H, i.e. [e~*^*, H] = Oo{£°°\t\), and it 
is associated with the spectral band aj-{z). The terminology was introduced 
in [Ne2]. Note, however, that Ranll is, in general, not an almost invariant 
subspace in the sense of [Ne2], since 11 need not have a limit as £ — > 0. 

Remark 2.2. Note that the growth condition on the gap in (IG)m is stronger 
than one would expect from the analysis in [NeSo] or [Tei]. Indeed, in both 
examples a gap which is bounded globally over phase space suffices to prove 
uniform adiabatic decoupling also in the presence of a Hamiltonian with 
principal symbol increasing linearly in momentum. More general, uniform 
adiabatic decoupling should hold whenever (IG)m is satisfied with 7 = m — p. 
Indeed (3) follows from the following proof with slight modifications under 
this weaker condition on the growth of the gap. However, this modified 
proof does not give vr G S^{e), a fact we will make use of in the following. 
To avoid further complications in the presentation, we decided to state only 
the stronger result for the stronger growth condition. 

Proof. We decompose the proof into two steps. 
Step I. Construction of the Moyal projector 

In general ttq is not a projector in the Moyal algebra, i.e. ttq # ttq ^ ttq. 
The following lemma shows that ttq can be corrected, order by order in 
e, so to obtain a true Moyal projector vr which Moyal commutes with H. 
Similar constructions appeared in the context of the Schrodinger equation 
several times in the literature [NeSo, BrNo, EmWe]. Our proof was strongly 
influenced by the one in [NeSo], but differs in relevant details, since we 
consider different symbol classes. It relies on the construction of the local 
Moyal resolvent of Hq{z). The construction of the global inverse of an elliptic 
symbol, often called the parametrix, is well known [DiSj , Fo, Ni] . 

Lemma 2.3. Assume either (IG)m or (CG). Then there exists a unique 
formal symbol 

j>0 
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such that '7Tq{z) is the spectral projector of Ho{z) corresponding to ar{z), with 
the following properties: 

(i) TT # TT = TT, 
(ii) TT* = TT, 

(iii) [H,7r]#:=H#7r-7r#H = 0. 

Proof. Wc give the proof under the assumption (IG)^. The proof under 
assumption (CG) is simpler, since all the symbols which appear belong to 

We first provide a constructive scheme for the special case where ar:{z) = 
{Ej-{z)} is an eigenvalue, which, at the same time, proves uniqueness of tt in 
the general case. It follows basically the construction as given in [EmWc]. 
The reason for including this scheme is that the aim of adiabatic perturbation 
theory is, in particular, to give an as simple as possible recipe for explicitly 
computing the relevant quantities. The inductive scheme for constructing 
TT in the special case cr^-^z) = {E^{z)} is much better suited for explicit 
computations than the general construction which will follow later on. 

Note that 7ro#7ro — ttq = 0{e) and [Hq^tto]^ = 0{e) and proceed by 
induction. Assume that we found tt^"^ = Yl]=o such that 

7r(") - 7r(") = e"+^ Gn+i + 0(e"+^) , (4) 

where, in particular, (4) defines Gn+i- Thus the next order term in the 
expansion 7r„+i must satisfy 

TTn+i TTo + TTq 7r„+i - 7r„+i = , 

which uniquely determines the diagonal part of tt^+i to be 

T^n+l = """"O Gn+l TTo + (1 - TTq) Gn+1 (1 " TTq) • (5) 

Since Gn+i = ttq Gn+i ttq + (1 — ttq) Gn+i (1 — T^o) follows from the fact that 
Gn+i is the principal symbol of e-"-i(^(") #^(") - ^W), cjW := tt^") + 
^"^^^n+i indeed satisfies (i) up to an error of order 0(6"+^). 

By induction assumption we also have that [iJ, 7r*^"^]# = 0{e"'~^^) and 
thus 

[H, J^'^y = s^+'Fn+i + 0(e"+2) . (6) 
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Hence, the diagonal part of vr^+i being fixed already, the off-diagonal part 
of TTn+i must satisfy [i7o,7rJ^+i] = —Fn+i- In particular, 

Ho{z) (7ro(z)7r„+i(2;)(l - 7ro(2:))) - {iro{z)Trn+i{z){l - 7ro(z))) -^0(2:) 

= -7ro(2)F„+i(z)(l-7ro(z)) (7) 



for all z G R . We first show that if (7) has a solution 7ro(-z)7r„+i(z)(l — 
7ro(z)) =: 7r^_P/(z), it is unique, i.e. that the kernel of the map Tr^^^{z) ^ 

[Hq{z) , Tr'f^^l^i (z)] restricted to Ran(l — 'Ko(z)) contains only zero. To see this 
let CTr(-2) := (sup(Tr(2;) — inf a"r(2;))/2 and note that, due to the gap condition, 
Hq{z) — a,:{z) is invertible on Ran(l — ^^{z)) with \\{Hq{z) — 'ar{z))^^{l — 
7ro(2;))|| < 2/diam((Tr(2;)). Hence 

[Ho{z),i,'^^hz)] = ^ [ifo(^) -^r(z),7rO^/(z)] = 
^ = {Ho{z)-a,{z))7:^^l{z){Ho{z)-a.{z)r' 

and therefore 

< \\iHo{z)-a.iz))Mz)\\ h^^i\z)\\ \\{Hoiz)-a,{z)r\l-7roim 
= C\\7rO^,\z)\\ 

with C < 1. Hence 7r^_Pi^(-z) = and we conclude that vr^+i is unique when 
it exists. 

In the special case that ar{z) = {Ej:(z)}, (7) can be solved, and one finds 

7ro7r„+i(l - TTo) = TTo Fn+i {Hq - E^)''^ (1 - TTo) . (8) 

Using that -Fn+i is the principal symbol of £~'^~^[H,uj^"'>]^, that ttq is the 
principal symbol of cu^"'^ and that o;^"-* satisfies (i) up to ©(e""*"^), one finds 
that 7roF„_|_i7ro = (1 — 7ro)F„4_i(l — ttq) = and thus that 7r("+-^^ defined 
through (5) and (8) satisfies (i) and (iii) up to ©(e"^^). 

Wc conclude that by induction we have uniqueness of vr in the general 
case, and an explicit construction for tt when ar{z) = {Ej.{z)}. The latter 
one involves four steps at each order: [a] Evaluation of Gn+i as in (4), [b] 
computation of as in (5), [c] evaluation of F^+i as in (6), [d] computa- 
tion of TT^+i as in (8) . 

We now turn to the construction of vr in the general case. Since the Moyal 
product is a local operation (it depends only on the pointwise value of the 
symbols and their derivatives) it suffices to construct tt locally in phase space 
and then uniqueness will liberate us from gluing the local results together. 
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Let us fix a point zq G M. . From the continuity of the map z ^ Hq{z) 
and the gap condition it foUows that there exists a neighborhood Uz^ of zq 
such that for every z G Uzq the set ar{z) can be enclosed in a positively- 
oriented complex circle T{zo) (independent of z) in such a way that T{zo) is 
symmetric with respect to the real axis, 

dist(r(zo), ct(^)) > ^Cg {pY for ah z G U^^ (9) 

and 

Radius(r(zo)) < C, sup {pY , (10) 
zeUzQ 

where Radius (r(;2o)) is the radius of the complex circle T = T{zq). The 
constant Cg in (9) is the same as in (2) and the existence of a constant 
independent of zo such that (10) is satisfied follows from assumption (G3). 
We keep a in the notation as a bookkeeping device, in order to distinguish 
the contributions related to the gap, although a = m. 

Let us choose any C ^ T and restrict all the following expressions to 
z G Uz^y There exist a formal symbol i?(C) - the local Moyal resolvent of H 
- such that 

R{0#{H-C,l) = l = {H-C,l)*R{Q onUz,. (11) 
The symbol i?(C) can be explicitly constructed. We abbreviate 

MO = {Ho-ar' 

where the inverse is understood in the ;B(7t!f )-sense and exists according to 
(9). By induction, suppose that i?(„)(C) = Z^"=o ^"'-^j(C) satisfies the first 
equality in (11) up to 0(£"'"'"^)-terms, i.e. 

i?(")(C) # (i? - CI) = 1 + £"+'^n+l(C) + 0(£"+2). 

By choosing Rn+i = -E^+i (Hq - Cl)~\ we obtain that = + 

e"''^^Rn+i satisfies the same equality up to C'(e"'"'"^)-terms. Then the formal 
symbol i?(C) = Sj>o^"'^7(C) satisfies the first equality in (11) which - by 
the associativity of the Moyal product implies the second one. 

Equation (11) implies that R{C) satisfies the resolvent equation 

R{0 - R{C') = (C - C') R{0 # RiC) on Uzo- (12) 



From the resolvent equation it follows - by using an argument similar to 
the standard one in operator theory [Kai] - that the symbol tt = X^j>o ^'''^j 
defined by 

•= i / ^) ^ e , (13) 
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is a Moyal projector such that [H,Tr]^ = on Uzq- Indeed, for every fixed 
z G and j € N, the map i— > -^^(C, z) is holomorphic in a neighborhood 
of the circle T{zo). Then r{zo) can be expanded to a shghtly larger circle V 
without changing the left hand side of (13) and we obtain 



I 

2^ 



where (12) has been used. The first equality in (14) follows by noticing that 
for every 7 G N^*^ 

(^) = ^ 92 R, {C,z)dC z€ U,, , 

and by expanding the Moyal product order by order in s. 

Since the circle T is symmetric with respect to the real axis one imme- 
diately concludes that vr* = vr, since R{C)* = R-iC) as a consequence of (11). 
From (13) it follows that vr Moyal-commutes with R{X) for any A G F. Then, 
by multiplying vr # -R(A) = -R(A) # vr by (i? — Al) on both sides, one obtains 
that iJ#7r = 7r#/i". 

Finally we have to show that ttj G Sp-"' for every j G N. From the Riesz 
formula (13) it follows that for every 7 G N^*^ one has 

||(a3'7r,-)(^)llB(7if) <27r Radius (F(zo)) sup ||(5Ji?,-) (C, ■ 

CGr(zo) 

According to (10) we are left to prove that 

(d^d^Rj) (C, z) < C^pj (j,)—^P-\^\P , a, /? G N^ i G N , 

(15) 

where Ca^j must not depend on zq. As for Rq, we notice that according to 
(9) one has 

and moreover, 

\\'^pRoiz)\\s(n() = l|-(-RoVpifo-Ro)(^;)||^(>^j) 



sup 

Cer(2o) 



< 



2 



Aj (p)— ||V,i7o(^)|lB(^,) 
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where the last bound follows from the fact that Hq G (recall that a = m). 
By induction one controls higher order derivatives and (15) follows for j = 0. 

Again by induction, assume that Rq, . . . ,Rn satisfy the bound (15). Then, 
by writing out 

£^n+l=(i2(n)(C)#(^-Cl)-l)„+i 

and using (82), one concludes that Rn+i = —En+iRo satisfies (15) with 
a = m. □ 



Step II. Quantization 

First of all, by resummation (Prop. A. 8) we obtain a semiclassical symbol 
TT : M^*^ X [0, eo) ^ S{Hi) whose asymptotic expansion is given by J2j>o ^'''^j- 
Then, by Weyl quantization, one gets a bounded operator tt G BiTi) (see 
Prop. A. 3) which is an almost-projector, in the sense that 

(i) 7f2 = 7f + 0_oo(£°°) 
(ii) TT * = TT 

(iii) [#,7r] = 0_oo(£°°) 



Notice that the assumption /> > is crucial in order to obtain (iii) for an 
unbounded H. 

In order to get a true projector we follow the idea of [NeSo] and notice 
that 11??^ — TT II = 0{£°°) and the spectral mapping theorem for self-adjoint 
operators imply that for each n G N there is a C„ < oo such that 

(7(7f) C [-C„e", C„£"] U [1 - C„£", 1 + (7„£"] =: (7^ U af . 

Hence one can define for e < l/(4Ci) 

27r 7|c-i|=i 

Then U'^ = U follows and we claim that 11 = tt + Oo{£°°). Indeed, 
9= [ XE{dX) = C'o(£") + / E{dX) = n + Oo(e") for all n G N , 
where £'(•) is the projection valued measure of ^. Finally notice that 

[#,n] = ^ / [^,(7f-C)-VC 

= -f / (?f-C)-M^,^](^-0-'rfC, 
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which impUes that 

\\[H,U]\\^^^^<C\\[H,n]\\^^^^=Oo{e^). 
This concludes the proof of the theorem. □ 

Essential self-adjointness of H. 

Since H is an hermitian symbol its Wcyl quantization H is symmetric on 
the invariant domain S{W^,B{Hi)) C 7Y. If if belongs to S^{e) then if is a 
bounded operator, and there is nothing to prove. 

In order to prove essential self-adjointness in the case H G Sl{e), we use 
an argument of [Ro]. The proof does not exploit the smallness of e and we 
therefore consider any e > 0. For s > let 

B±siq,p) = iHoiq,p)±isl)-\ 

which, according to Proposition A. 5, belongs to Si{B{H{)). Moreover 

{H ± isl) # = 1 + eS±s , 

where S±s G S^{£), since H G S'J(e) and B±s G S^{e). After Weyl quantiza- 
tion we obtain that 

^ ^ ^ ^ C 

{H ± isl) B±s = 1 + £S±s with H-SisHg^^) < j^, 

the latter bound following (for s large enough) from Proposition A. 3 and 

from estimating the Frechct semi- norms of S±s. Essential self-adjointness of 
H on the domain <S follows, if we can show that Ker(i/* ibis) = {0} for some 
s > 0. For this let if G Ker(if* it is) and ip ^ S. Using B±S C S, we obtain 

= {{H* ^ is)^, B±i^) = {if, {H ± is)B±i^) = (1 + sS±s)^) . 

Since ||e5'±s|| < 1 for s large enough, (1 -|- £S±s)S is dense in H and hence 
if = follows. 

3 Reference subspace and intertwining unitaries 

The fact that the subspace associated with an isolated energy band decouples 
from its orthogonal complement up to small errors in e leads immediately to 
the following question. Is there a natural way to describe the dynamics of the 
system inside the almost invariant subspace Ranll? The main obstruction 
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for such a simple description is the fact that the subspace RanlT depends on 
e and is not easily accessible. Even worse, in general the limit lim^^o n does 
not exist, meaning that Ranll is not even close to an e-independent subspace. 
In order to obtain a useful description of the effective intraband dynamics we 
thus need to map Ranll to an easily accessible and e-independent reference 
subspace. 

^Prom the continuity of z ^ Hq{z) and the gap condition it follows that 
there is a subspace /Cf C 7i{ independent of {q,p) such that the subspaces 
Ran7ro(g,p) are all isomorphic to /Cf. Let TTr be the projection on /Cf, then 
Hi. := 1(8) TTr (= TTr) will serve as the projector on the reference subspace JC := 
Ranllr. Of course /Cf is highly non- unique and a convenient choice must be 
made in concrete applications. 

Once the reference Hilbcrt space is fixed we next chose a unitary operator 
valued smooth function uq{z) which pointwise in phase space intertwines 
7ro{z) and VTr, i.e. 

Uo{z) Tro{z) Uo{z)* = TTr . (17) 

The existence of such a smooth map follows from a bundle-theoretic argu- 
ment given at the end of this section. Again uo{z) is not unique and must 
be chosen conveniently. We will see in Section 6 that there is an optimal 
choice for uo{z), which reflects the physics of the problem. 

We cannot prove that it is possible to choose uq in S^{B{Tlf)). Indeed, 
relation (17) does not imply any bound at infinity on the derivatives of uq, 
as can be seen by multiplying uq with a highly oscillating phase. Hence we 
assume that uq is in Sp{B{'Hi)), as will be the case in the physical examples. 

In the following U{H) will denote the group of unitary operators over H. 

Theorem 3.1. Assume either (IG)m or (CO) and that there exists aU{Ti.f)- 
valued map uq G Sp(B{7if)) which satisfies (17). Then there exist a unitary 
operator U € B{H) such that 

UUU* = U, (18) 
and U = u + Oo{£°^), where u x Y^j>Q£''uj in S^p{e) with principal symbol 

Uq. 

Remark 3.2. In [NeSo] the Nagy transformation (22) is used in order to 
map Ranll to the £-independent subspace RanTTQ. This is possible because 
in their application the symbol ttq depends only on q and, as a consequence, 
ttq is a projector satisfying ||n — vfoH = 0{e). However, in general ttq depends 
on q and p, see Section 6 and [PST2] for relevant examples, and the mapping 
to the reference space becomes more subtle. 
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Proof. Step I. Construction of the Moyal unitaries. 

Again uq fails to be a Moyal unitary (i.e. Uq^uq ^ 1) and to intertwine tt 
and TTr- However, the following lemma shows that uq can be corrected order 
by order to reach this goal. The idea of constructing a pseudodifferential 
operator which is almost unitary and diagonalizes a given pseudodifferential 
operator has a long tradition, cf. [Ni] Section 7 and references therein, and 
was applied in different settings many times, e.g. [Ta, HeSj]. 

Lemma 3.3. Assume either (IG)m or (CG) and that there exists aU{TC{)- 
valued map uq G S^p{B{l-if)) which satisfies (17). Then there is a formal 

symbol u = '}2j>o^''''^j > % ^ Sp-"^{B{H{)), such that 

(i) u*#u = l and u#u* = l, 

(ii) U # TT # U* = TTr , 

where ir is the Moyal projector constructed in Lemma 2. 3. 

Remark 3.4. We emphasize that - as opposed to the Moyal projector tt 

appearing in Lemma 2.3 - the Moyal unitary u is highly non-unique even 
for fixed uq. As it will follow from the proof, all the possible choices of 
Moyal unitaries intertwining tt and VTr with prescribed principal symbol uq 
are parametrized by the antihermitian Moyal symbols which are diagonal in 
the TTr-splitting. 



Proof of Lemma 3.3. Observe that uq satisfies (i) and (ii) on the principal 
symbol level. We proceed by induction and assume that we found u^"'^ = 
^^=(j£-^Uj satisfying (i) and (ii) up to Oie"^^^). We will construct lin+i such 
that «("+^) = +£"+i«„+i satisfies (i) and (ii) up to 0{e"+'^). To this 
end we write without restriction 

Un+l =■ (On+l + bn+l) Uq , 

with ttn+i hermitian and 6^+1 anti-hermitian. By induction assumption we 
have 

nH#nW*-l = £"+iA„+i + 0(£"+2) 
Thus Un+i has to solve 

UoU*i+l + Un+lU^ = -An+1, ^^^^ 
■Wj^+l Uo + Uq Un+l = -An+1 . 
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The first equation in (19) fixes a„+i = since An+i is hermitian as 

it is the principal symbol of e^^~^{u^"'^ ^^n^")* — 1). The second equation 
in (19) is then also satisfied, since the compatibility equation uq An+i = 
An+i Uq follows from 

# («(-)*#«(") - 1) = -i^(nW - 1) 

by noticing that uq ^n+i (resp. ^n+i uq) is the principal symbol of the l.h.s 
(resp. r.h.s). 

Note that (19) puts no constraint on 6n+i Siiid we are left to determine 
it using (ii). Let w^'^^ = u^'^^ + e""*"^ a„+i uq, then by induction assumption 

#7r - TTr = + 0(£"+2) 

and thus 

#7r#n("+l)* - TTr = (S„+i + [6„+l,7rr]) + 0(£"+') . 

Hence we need to find an anti-hermitian bn+i satisfying 

Bn+l + [bn+l,'^r] =0, 

which is given by 

6„+l = [vTr, Bn+l] , (20) 

provided that -Bn+i is hermitian and off-diagonal in the TTr-splitting, i.e. 
TTr Bn+l TTr and (1 — TTr) Bn+l (1 — TTr) vanish. This follows by noticing that 
Bn+i is the principal symbol of e~("'+^) (■«;("■) # tt # — TTr) and then 

(1 - 7r,)i?„+l (1 - Vr,) ^(1 - TTr) (^«(") #7r#t/;(")* - TTr) (1 - TT,) 

= - TTr) (^«(") # TT # ^(")*) (1 - TTr) 

= ^ (^'^"-''^i^n+l (^^"^ # ^# ^^"^*) + 0(e"+2)) 0, 

where for the last equality we inserted 1 — tTi- = w^'^^ 7^ (1 — vr) ^ -|- 
e^^+^Sn+i + 0(6^^+2) and used that solves (i) up to 0{e''+^) and that 
TT is a Moyal projector. A similar argument shows that tt^ Bn+i tTj. vanishes 
too. Note also that (20) fixes only the off-diagonal part of bn+i and one 
is free to choose the diagonal part of arbitrarily, which is exactly the 
non-uniqueness mentioned in Remark 3.4. 

It remains to show that the assumption uq G implies that Uj belongs 
to Sp-"^. Assume by induction that u^"'^ G Mp(£). Then the formula 

an+l = —An+l = -I - 1) 

Z Z \ / n+1 
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shows that a„+i belongs to Sp ^"''^^^^ as it is the (n + l)-th term of an el- 
ement of Mp{e). By Proposition A. 4, a„+i G S^*-""*"^^^ as well. Anal- 
ogously we have that Bn+i € Sp '•"^^^'^ by induction assumption, there- 
fore bn+i € Sp and thus bn+i uq G Sp ^"^^^^^^ which finally gives 



Step II. Quantization 

Now let u denote a resummation of the formal power series u = Ylj>o ^'''^j 
in S'^(e) (sec Prop. A. 8). Then, by Wcyl quantization, one gets a bounded 
operator u G B{7i) (see Prop. A. 3) such that: 

(i) = 1 + 0_oo(e°°) and uu* = 1 0_oo(e~) 

(ii) MvfS* = n, -t-0_oo(£°°). 

As a first step we modify u by an C'o(e°°)-term in order to get a true unitary 
operator U G U{J-L) (which, in general, does not correspond to the Weyl 
quantization of any semiclassical symbol). Let 

U = u{u*u)~^ . (21) 

Notice that u*u is a self-adjoint positive operator which is Oo(e°°)-close 
to the identity operator. Then {u*u)~'2 is well-defined and again Oo{e°°)- 
closc to the identity operator. Hence (21) defines a unitary operator which 
moreover is C'o(e°°)-close to u. 

Finally we modify U in order to obtain a unitary which exactly inter- 
twines Ilr and n. Since ||C/ IT C/* — Ilrll < 1 for e sufficiently small, the Nagy 
formula as used in [NeSo] 



W := 



1 - (uuu* -U, 



UUU*U^ + {1-UUU*){1-U^) (22) 



defines a unitary operator W G U(H) such that W UUU* W* = Uj. and 
W = 1 + Oo(e°°)- Thus by defining U = W U one obtains (18), with the 
desired properties. □ 

Remark 3.5. We sketch how to prove the existence of a smooth map uq 
satisfying (17). Given 

£; = |(z, V') G M^'^ X Hf : V e Ran7ro(z)| 
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the map He E ^ R {z, tl:) ^ z defines a fibration of Hilbert spaces over 
the base space M^*^. 

The fibration is locally trivial. Indeed for any zq G M^*^ there exists a 
neighborhood Uz^ such that ||7ro(2;) — 7ro(-2;o)|| < 1 for any z G Uzq-, so that 
the Nagy formula 

_ 1 

w{z) = 1 - {itq{z) - 7ro(zo))^ ^ [7ro(^;)7ro(^;o) + (1 - t^o{z)){^ - T^oizo))] 

locally defines a unitary operator w{z) such that w{z)* ttq{ z)w{z) = 7ro(2;o)- 
A local trivialization of the fibration is then explicitly given by 

e : Y\.'^^{Uzq) Uzo X Ran7r(zo) ^ l^zo x ACf 

i-^' {z,w{z)'iIj) hh- {z,(j){zo)w{z)'(p) 

where we use the fact that there exists a unitary operator 4>{zo) : Ran7r(2:o) — > 
/Cf. The existence of ^(^o) follows from the fact that the dimension of 
Ran7r(2;o) is independent of zq, but the map zq h- > ^(2:0) may be a priori 
even discontinuous. 

Moreover one can check that any two such trivializations are IA{1C^)- 
compatible, and the previous data define a linear W(AIf)-bundle. 

Since the base space is contractible, the bundle is trivial and the asso- 
ciated principal ^(/Cf) -bundle (i.e. the bundle of the orthonormal frames) 
admits a global smooth section. This implies the existence of a smooth map 
uo : M^*^ U{ni) such that (17) holds true. 



4 Adiabatic perturbation theory 
4.1 The effective Hamiltonian 

In the previous section we constructed a unitary U on H which exactly 
intertwines the almost invariant subspace Ranll and the reference subspace 
/C = Ranllr. U and 11 are C'o(e'^)-close to pseudodifferential operators with 
symbols u and tt both in Sp{e). 

We define the effective Hamiltonian h as the quantization of a resumma- 
tion h of the formal symbol 



h = u#H#u*. 



(23) 
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Recall that wc do not distinguish semiclassical symbols and formal symbols 
in the notation. The following theorem is the basis for the adiabatic per- 
turbation theory, as it relates the unitary time-evolution generated by the 
original Hamiltonian H to the one generated by the effective Hamiltonian h. 

Theorem 4.1. Under the assumptions of Theorem 3.1, one has that h G 
S'^ie) and h is essentially self-adjoint on S. Furthermore 

[/i, nr]=0, (24) 

e-'^^ -u*e-'^^u = Oo{e'^\t\) (25) 

and 

^-iHt _ ^* ^-iht ^ ^ Oo{s°°{l + \t\)) . (26) 



Proof. Since u G Sp{£) and H G S^{e), the composition rule for semiclas- 
sical operators (see Prop. A. 6) yields h G S^{e) and thus hj G S^~-"^. 

Let h ■.= uHu*. Since u* is bounded with bounded inverse, one finds, by 
checking definitions, that h is self-adjoint on u*~^D{H) and that h is essen- 
tially self-adjoint on u*~^S. According to Equation (8.10) in [DiSj], which 
generalizes to 6('Hf)-valued symbols, u*^^ G OPS^{e) and thus u*^^S = S. 
Hence <S is a core for h and, since h — h & BiTi), the same conclusions hold 
for h. 

Next observe that, by construction, [/ij,7rr] = for all j G N and thus 
[/ij,7rr]^ = because TTr does not depend on {q,p) G M^*^. Hence [/ij,nr] = 
and thus (24) follows. 

For (25) observe that 

e-'"*-u* e-'^* u = -i e-'"* ds e'^' (^Hu*-u*h^ e'"''' u = Oo{e^\t\) , 

since, by construction, {Hu* — u*h) = 0_oo(£°°). Finally (26) follows from 
(25) using C/-u = Oo(e°°). □ 

Remark 4.2. It might seem more natural to define the effective Hamiltonian 

as 

iJeff = UUHUU* -\-U{l-U)H{l-U) U* . 

Clearly one should have H^s — h = 0{e°°) in some sense. However, if 
H is unbounded, this closeness does not follow in the norm of bounded 
operators from our results, since U need not be a semiclassical operator. 
As a consequence no asymptotic expansion of i?eflf in the norm of bounded 
operators would be available. 
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In the remainder of this section we will study the finite order asymptotic 
approximations 

n 
j=0 

to the effective Hamiltonian h. By virtue of (24), we can, whenever ap- 
propriate, restrict our attention to the reduced Hilbert space /C = Ranllr. 
Furthermore we define n^") = J2^j=o uj and obtain a finite order expansion 
of the unitary U as \\U - ||bc^) = C'(£'*+^). 

Our main interest are approximations to the solution of the time-depen- 
dent Schrodinger equation 

■ dipt £> , 
^^ = ^'^* 

over times of order e~^T, where r does not depend on s and G N is 
arbitrary. Starting with (25) on the almost invariant subspace we obtain 

e-'^^U = u*e-''^^'U,u + Oo{e°°\t\) 

= n^")* e-^*'^^"^'^ nru(") + (1 + |t|)Oo(£''+^) , \t\ < £~V(^7) 

where p{n -|- + 1) > m is assumed in order to have h — ^("■+*^) G B{H). 
Hence, given the level of precision e" and the time scale e"'^, the expansion 
of h must be computed up to order hn+k and the expansion of U up to order 
Un- Put differently, in order to improve the error, a better approximation 
to the unitary transformation is necessary. On the other hand, in order to 
enlarge the time-scale of validity for the space-adiabatic approximation, only 
the effective Hamiltonian h must be computed to higher orders. 

Specializing (27) to n = and k = 1, one obtains the leading order 
solution of the Schrodinger equation as 

e-'^*U = u*o e-^('^o+^'^i)* Ur m + {I + |r|)Oo(e) , \t\ < e"V , (28) 

where m < 2p. Here the choice of = 1 corresponds to the macroscopic 
or semiclassical time-scale t/e. On this time-scale the effective dynamics 
g-ift.t/£-Q^ on the reference subspace is expected to have a nice semiclassical 
limit, under suitable conditions on h. 

Note that one can replace in (27) and analogously in (28) r by e~^T and 
obtains 

e-'"^U = e-*^'"+'' H, S^*^) + {1 + |t|)Oo(£"+^-^) , \t\ < £-(*^+^V . 

(29) 
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Thus one can enlarge the time-span for which the approximation holds with- 
out the need to compute further terms in the expansion. The price to be 
paid is a larger error, of course. 

We emphasize that (27) and (28) are purely space-adiabatic expansions 
with no semidassical approximation invoked yet. As a consequence one ob- 
tains uniform results and a simple bound on the growth of the error with 
time. Note in particular that the space-adiabatic approximation holds on 
time-scales far beyond the Ehrenfest time-scale, the maximal time-scale for 
which semidassical approximations are expected to hold. For some partic- 
ular cases semidassical expansions of the full propagator e~*^*/^ have been 
derived directly, e.g. in the context of the Dirac equation [Ya, BoKe2]. These 
expansions hold, in general, only for short times, in the sense that they must 
be modified each time a caustic in the corresponding classical flow is en- 
countered. More important, the clear separation of the space-adiabatic and 
the semidassical expansion is not maintained, which is a severe drawback, 
since in many physical situations the space-adiabatic approximation is valid 
to high accuracy, while the semidassical approximation is not, cf. Section 6. 
On the other hand, a semidassical expansion of the right hand side of (28) is 
straightforward in many interesting cases, as will be discussed in Section 5. 

In parentheses we remark that the space-adiabatic approximation can be 
used also in the time-independent setting, i.e. to estimate spectral properties 
of H. If one is able to compute eigenvalues of h^"'^ up to errors of order o(£"), 



If, in addition, one knows from sonic a priori arguments that H has pure 



explicit knowledge of U is needed as long as the interest is in approximate 
eigenvalues only. For example, the scheme just described can be applied to 
the time-independent Born-Oppenheimer theory, where one is interested in 
the low lying spectrum of a molecule. The standard approaches to the time- 
independent Born-Oppcnheimer approximation [CDS, Hai, KMSW] yield in 
some respects mathematically stronger results. However, our scheme suffices 
for estimating asymptotic expansions of eigenvalues and is simpler to handle, 
in general. 



it follows that 
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4.2 Leading order terms in the expansion of the effective 
Hamiltonian 

We turn to the explicit determination of the leading order terms hj in the 
expansion of h using (23). Of course, in concrete applications only H and 
uq are given explicitly, while the higher order terms in the expansion of u 
must be calculated using the construction from Section 3. For a general 
Hamiltonian H such a program is feasible only for the terms ho, hi and 
possibly /i2, which will be our concern in the following. 

The principal symbol of h is given by 

ho = Uq Ho Uo . 

Higher order terms can be obtained using (23). The double Moyal product 
becomes rather awkward to handle, and alternatively we proceed inductively 
by observing that 

u#H -ho#u = ehi#u + 0{e'^) = ehmo + 0{e'^) , (30) 

with the subprincipal symbol on the left hand side being 

{u#H - ho#u)i = uiHq + uqHi - houi + {uo # -ffo)i - {K # uo)i ■ (31) 

Recall the notation a^^b = Yl'^o # expansion of the Moyal 

product, see the Appendix. Combining (30) and (31) one obtains 

hi = [uiHo + uoHi - houi + {uq # i?o)i - (^o # ""0)1 ) ^^o • (32) 

The expression (32) further simplifies if one specializes to the case where 
<^r{q,p) = {Er{q,p)} consists of a single eigenvalue of Ho{q,p) and one 
projects on the relevant subspace, 

TT^hiTTr = TTr{uo Hi Uq + (tlQ # i^o)l ^^0 " i^* # ""0)1 Wo)^r • (33) 

The right hand side has the nice property to be independent of ui and thus 
to depend only on known quantities. 

Along the same lines and under the same condition on ar{q,p), one com- 
putes 

7rr^27rr = Trr(^uoH2 + uiHi — hiUi (34) 
+ {ui # Ho)i + {uo # Hi)i - {E^ # ui)i - {hi # uo)i 

+ {uo # Ho)2 - {E* # '^0)2) ullTr . 



168 



Space- Adiabatic Perturbation Theory 



Again, (34) does not depend on U2 for the special case under consideration, 
but it does depend on ui, which must now be computed using the construc- 
tion from Section 3. 

Although (34) looks still rather innocent, in general, it requires some 
work to compute it explicitly. This is partly because the second order expan- 
sion of the Moyal product in (34) tends to become rather tedious to obtain. 
But, in general, also the determination of ui is nontrivial. To convince the 
reader, we state without details that the construction from Sections 2 and 3 
yields 

«i = ^o( - ^{uo,Uo} + [uoT^i^ Uq,t^t] +^[{{uo,no}Uo + Uo{Tro,UQ}),7rr]^ , 

(35) 

with 

7r°° := 7ro7ri(l - ttq) + (1 - 7ro)7ri7ro , 

where we used that (a#6)i = — |{a, 6}. Recall the definition (83) of the 
Poisson bracket {•,•}. 

To compute tti from the given quantities one has to use the construction 
explained in Section 2. One finds 

7r?° = '-{RoiE.)il-TTo){Ho + E„TTo}Tro 
+ vro {vro, Hq + E,.} Ro{E^) (1 - ttq)) 
+ 7roi/ii?o(^*)(l - vro) + i?o(^*)(l - vro)i^ivro , 

where Ro(Ey-){l — ttq) = {Hq — £^r)~^(l — ttq) is uniformly bounded because 
of the gap condition. For sake of completeness we mention that tti = ttP^ -|- 
§{7ro,7ro} in this case. 

For the higher orders in the expansion of h we only remark that, in 
general, /t„ depends on and h^^^^\ In the special, but interesting 

case of an isolated eigenvalue Ei-{q,p), hn depends only on u^^~^\ H^^^ and 
/i^""^) and is thus considerably easier to obtain. 

Remark 4.3. Note that in the case of a;-{q,p) = {E,.{q,p)}, not only the 
principal symbol ho{q,p) = Ej:{q,p)l-}if, but also the subprincipal symbol 
hi{q,p) as given by (33) is well defined regardless of the gap condition, 
provided that the spectral projection TTo{q,p) is sufficiently regular. Indeed, 
it can be shown, at least in some special cases, that there is still adiabatic 
decoupling to leading order and an effective dynamics generated hy ho + e hi 
without a gap condition [Tei], [Te2]. 

To get even more explicit formulas for hi and /12, note that in most 
applications one has no naturally given transformation uq. Instead one 
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chooses a suitable basis {'>pa{QiP)}ael of Ran7ro((j',p) and defines uo{q,p) = 
X^ae/ \Xa){i^aiQ,p) \ + f{<liP), where the vectors Xa form a basis for RanTTr 
and r{q,p) is some arbitrary unitary intertwining Kamro{q,p)-^ and Ran7r;f . 
TTj. hj{q,p) TTj. is independent of the choice of the unitary r{q,p) for all j G N. 

We remark that such a basis {ipaiQ,p)}aei of global smooth sections of 
the bundle over M^'^ defined by TTo{q,p) always exists, since M^'^ is contractible 
(see Remark 3.5). However, we are not aware of a proof which insures 
uq € Sp . The situation changes completely, once one considers local domains 
in the base space which are not contractible. Then it might become necessary 
to chose as reference space the space of sections of a globally nontrivial 
bundle. 

Assuming that (Tr{q,p) = {Ej-{q,p)} consists of a single eigenvalue of 
Ho{q,p) of multiplicity £ (including £ = oo), we obtain the £ x ^-matrix 

TTr /i^^-*(Q',p) TTr aS 

= ' ^^'^ Xp)=Er Sa/3 +ehiaP, (36) 

with 

= (V'a,^^lV'/3) -i(V'a,K,V'/3}) " ^ (Va, {(^0 - ^r), V'/s}) • (37) 

The indices a and /3 are matrix-indices, both running from 1 to £. Equations 
(36) and (37) are one of our central results. They are still of a simple form 
and mostly suffice to compute the basic physics. The first term in (36) 
is referred to as Peierls substitution and the first order correction carries 
information on the intraband spinor evolution. E.g., as will be discussed 
in Section 6, for the Dirac equation hi governs the spin precession. The 
reason for the particular splitting of the terms in (37) will be discussed in 
Section 5. Here we only remark that the second term in (37) is related to 
a "generalized" Berry connection. We omit the analogous formula for h2a/3, 
since it is too complicated to be helpful. 

4.3 Born-Oppenheimer type Hamiltonians 

An instructive example to which formula (37) applies are Born-Oppenheimer 
type Hamiltonians of the form 



HBo{q,p) = ^p^lHf+Viq), 



(38) 
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V G S^{B{Tl{)), with an electronic energy band ej-{q) of constant multiplicity 
£, i.e. V{q)Tro{q) = er(g)7ro((/). Adiabatic decoupling for Born-Oppenheimer 
type Hamiltonians is established with exponentially small errors by Martinez 
and Sordoni [MaSo] , see also [So] . Their result partly triggered our interest to 
develop a general theory. Exponentially accurate coherent state solutions for 
Born-Oppenheimer type Hamiltonians have been constructed by Hagedorn 
and Joye in [HaJo2]. 

Note that the quadratic growth of H^Q{q,p) as a function of p prevents 
from applying the general results directly. As to be discussed in Section 4.5, 
energy cutoffs need to be introduced. For the moment we ignore this problem 
and proceed by working out the perturbative scheme formally. 

We fix arbitrarily an orthonormal basis {tpaiQ)}a=i of Ran7ro(q') depend- 
ing smoothly on q which then satisfies H-Qo{q,p)ipa{q) = Er{q,]:))il!a{q) with 
Er{q,p) = ^P^ + fir (9) for 1 < a < ^. Only the second term of our formula 
(37) contributes and yields 

hlaf}{q,P) = -ip- {ll^a{q),yqlkp{q)) =■ -P-Aaf}{q), 

which is well known in the case of a nondegenerate eigenvalue, [ShWi, LiWe, 
TeSp]. Aap{q) has the geometrical meaning of a gauge potential, i.e. coef- 
ficients of a connection on the trivial bundle x C^, the so called Berry 
connection. As mentioned already, a more detailed discussion of the origin 
of the Berry connection will be given in Section 5. 

For the Born-Oppenheimer Hamiltonian the calculation of /t2a/3 is still 
feasible without much effort and the result is 

1 ^ 1 

-{p-Vq^a,Ro{Er)p-Vq^f3). (39) 

Recall the definition of Ro{Ey.) = {Hq — Ej.)~^{l — ttq), which reduces to 
RQ{Ey.){q) = {V{q) — er(g))~^(l — 7ro(g)) in the present case. Although we 
omit the details of the computation leading to (39), we shortly describe how 
(34) relates to (39). Since Hi = and H2 = the corresponding terms in 
(34) do not contribute. Since uq and ttq are functions of q only, the second 
term in (35) is the only one contributing to ui, and thus the third term in 
(34) also vanishes after projecting with the tTi's from outside the brackets. 
The last two terms in (34) cancel each other. The seventh term in (34) yields 
the first term in (39) and the fourth and sixth term in (34) combine to the 
second and third term in (39). In particular the calculation yields for the 
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symbol of the unitary 

e 

u*Bo(.q,p)7^. = Yl {\MQ))+i^MEr){q) \p ■ VqMQ))) {Xa\ + 0{e^) . 

a=l 

Thus the symbol of the second order effective Born-Oppenheimer Hamilto- 
nian reads 

1 / \ 2 

+y (V,Va(g), (1 - Mq)) ■ v.VM^)) (40) 
-e^{p- VgMq), MErM p ■ VqMi)) + <^(^^) > 

where the first term from (39) nicely completes the square to the first term 
in (40). Note that the third term on the right side of (40) depends on q only 
and was interpreted in [ShWi] as a geometric electric potential in analogy to 
the geometric vector potential A{q). 

In the special case of a nondegenerate eigenvalue Cr and a matrix- valued 
Hamiltonian H, (40) reduces to the expression obtained by Littlejohn and 
Wcigcrt [LiWe]. They also remark that the previous studies [ShWi, AhSt] 
of the expansion of the effective Born-Oppenheimer Hamiltonian missed the 
last term in (40). This strengthens our point of the usefulness of a general 
and systematic space-adiabatic perturbation theory. 

The full power of our scheme is in force in cases where RauTTo is de- 
generate and depends both on q and p, since then the known techniques 
[LiFl, LiWe, NeSo, MaSo] cannot be applied. The simplest example of this 
kind is the one-particle Dirac equation with slowly varying electric and mag- 
netic potentials, which will be discussed in Section 6. 



4.4 The time-adiabatic theory revisited 

With little additional effort our scheme can be applied even to the time- 
adiabatic setup. As for notation, wc replace the phase space x by 
Mt X in the following. Given a Hilbert space Ti and family H^{t), t € M 
of self-adjoint operators such that H^{t) =: H(t,r],e) G S^{£,B{H)), the 
solutions of the equations 

isdtU'it, s) = H%t)U'{t, s), s G M , (41) 

define a unitary propagator. A unitary propagator is a unitary operator- 
valued map U{t, s) strongly continuous in t and s jointly, such that 

Uit, t) = In and Uit, r)U{r, s) = U{t, s) 
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for any r,s,t G M. In particular we have that i7^(t, 0)^o solves the time- 
dependent Schrodinger equation 

ie^^m=H'{t)m- (42) 

for any ipo ETC. 

It is assumed in addition that Ho{t), the principal symbol of H'^{t), has 
a relevant part (j^{t) of its spectrum, which is separated by a gap from the 
remainder uniformly for t G M. As before we denote the spectral projection 
on (7r(t) by 7ro(t). 

The following theorem is a variant of the time- adiabatic theorem of quan- 
tum mechanics [Ka2, ASY, JoPf, Nei], however formulated in the language 
of adiabatic perturbation theory. Sjostrand first recognized the usefulness 
of pseudodifferential calculus in this context [Sj] and we are grateful to G. 
Nenciu for pointing this out to us. We remark that the proof below can 
be adapted to the case of a time-dependent operator-valued classical sym- 
bol H{q,p,t), as - for example - the Dirac Hamiltonian or the Pauli-Fierz 
Hamiltonian with slowly varying time-dependent external potentials. 

Theorem 4.4 (Time-adiabatic theorem) . LetH{t) andar{t) be as above. 

(i) Decoupled subspace. There exists a family of orthogonal projectors 
U{t) such that n(-) G S°{e,B{n)), U{t) - 7ro(t) = Oo{e) and 

U(t, s)* U{t) U{t, s) = U{s) + Oo{e°^\t - s\) (43) 

uniformly for G M. Whenever dfH{t) = for some t G M and all 
a G N, then U{t) = 7ro(t). 

(ii) Intertwining unitaries. There exists a family of unitaries uo{-) G 
C{j°(M, BiTi)) with uo{t) 7ro(t) u^it) = 7ro(0) =: TTr and a family of uni- 
taries U{-) G S^{e,B{n)) such that 

U{t) U{t)U*{t) = TTr and U{t) - uo{t) = Oo{s) . 

(iii) Effective dynamics. There exists a family of self-adjoint operators h(t), 
hi-) G S^{e,B{n)), such that 

[h{t),7rr]=0 for all t eR (44) 

and the solution of the initial value problem 



isdtUeS{t,s) =h{t)Ues{t,s), SGR, Uesit,t) = lH 
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satisfies 

U{t,s) = U*{t)U^s{t,s)U{s) + Oo{e°°\t - s\) . (45) 
The asymptotic expansion of h{t) in B{7i) reads 

oo / 

Kt) >^ E^" E ^jit)H,{t)uUt) (46) 

n=0 \j+k+l=n 

+ ^ E {uj{t)ul{t)-uj{t)uUt)) 

j+k+l=n 

where £"'-ffn(i) is the asymptotic expansion of H{t) in B(l-L) and 
^n£^Un is the asymptotic expansion oflA{t) inB{H). 

Before we turn to the proof we remark that, for ar{t) = {er{t)} and 
{fa{t)}a=i orthonormal basis of Ran7ro(t), the effective Hamiltonian in- 
cluding second order reads 

where -Ro(er) = {H{t) — er(t))^^ (1 — 7ro(t)). For the unitary U{t) one finds 

e 

U*{t)n, = E {\^a{t)) + ieRo{er){t) \^a{t))) {'fam + O(e') • 

a=l 

Proof. In order to apply the general scheme developed in the previous 
sections it is convenient ~ in analogy with the extended configuration space 
in classical mechanics - to introduce the extended space K, = L^{M.,H) = 
H. dt and to define the extended Hamiltonian 

K = -iedt + H{t) 

which is self-adjoint on the domain T>{K) = i7^(R,7^) C /C. By following 
Howland [Ho], we notice that the unitary group e"*'^'^, a G M, is related to 
the unitary propagator (41) through 

(e-'^^if) (t) = U{t, t - a)tp{t - a) . (47) 

Moreover, the unitary group e~^^" can now be studied by means of the 
techniques developed in the previous sections, since K is nothing but the 




174 



Space- Adiabatic Perturbation Theory 



Wcyl quantization of the operator- valued function K(t,r)) = r) + H(t), and 
K belongs to 

By assumption K e satisfies assumption {Ga.p)a with a = 0. How- 
ever, because of the simple dependence of K(t, rj) on ry, the conclusion of 
Theorem 2.1 and 3.1 hold still true in a sense to be made precise. 

Indeed, by following the proof of Lemma 2.3 one obtains a semiclassical 
symbol vr G SQ{e,B{H)), depending on t only, such that [K,Tr]^ x in 
Sq^e). On the other hand, 

[K, vr]^ = [H, tt]^ + [r], n]^ = [H, vr]^ - ie {dtrr) 

where the last equality follows from the fact that [r/, 7r]_^ is the symbol of 
[—iedt,7r{t)] = —ie {dtir) {t). Since both [H, vr]^ and dti^ belong to iSg (e), one 
concludes that the asymptotic expansion [i^, vr]^ x holds true in Sq{£), 
and hence [i^, tt] = Co(£°°)- Finally one defines 

m = i / - cr'dc 

27r 7|c-i|=i 

and finds n(-) e 5°(e,B(7^)), U{t) - 7r(t) = Oo{e°°) and [e-^^^,n] = 
Oo{e°°\a\) as in Section 2. Together with (47) this implies 

ess sup \\U {t,t- a)* U{t)U{t,t- a) -U{t-a)\\,g(^. = Oo(e~kl)- 

However, since n(t) and U(t,s) are continuous functions of t, the pointwise 
statement (43) follows. 

For no(t) one can use for example Kato's construction [Ka2] and define 
uo{t) as the solution of the initial value problem 

j^u*o{t) = [Mt),Mt)]u*o{t), ^a(o) = i- 

Clearly uo{t) belongs to S^{B{'H)). Notice that the same construction does 
not work in the multidimensional case, since the evolutions in different di- 
rections do not commute. U can be obtained as in Section 3, where the 
fact that Tr{t) and uo{t) both depend on t only and not on rj simplifies the 
construction considerably and yields, in particular, a fibered unitary U{t). 

As in the general setting let the effective Hamiltonian be defined as a 
resummation of 



kirj, t,e) = {u#K# u*) (r?, t, e) =: r? + h{t, e) , 
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with the exphcit expansion (46). According to Theorem 4.1 we then have 

which imphes according to (47) that 

esssup\\U{t,t-a)-U*{t)Uesit,t-a)U{t-a)\U^. = Oo{e°°\a\) . 
teM ^ ' 

The pointwise statement (45) follows again from the continuous dependence 
on t of all involved expressions. □ 

4.5 Energy cutoff 

The Born-Oppenheimer type Hamiltonians as well as many other physically 
relevant Hamiltonians do not satisfy the general assumptions we imposed in 
Sections 2 and 3. This is so for two reasons. First of all they are quantizations 
of symbols taking values in the unbounded operators. Secondly, the gap does 
not increase as fast as the Hamiltonian for large momenta, e.g. quadratically 
in the Born-Oppenheimer setting. The first problem is purely technical and 
the domain questions which arise have to be dealt with case by case. The 
second problem causes a qualitative change in the sense that the adiabatic 
decoupling is no longer uniform, as can be seen from the construction of the 
almost invariant subspace in Section 2. To deal with the second problem 
one therefore needs a cutoff for large momenta. There are basically two 
ways to implement such a cutoff. One possibility is to directly cut off large 
momenta as was done in [TeSp, SpTe], but then one needs to control the 
times for which no momenta exceeding the cutoff are produced under the 
dynamics. However, for a large class of Hamiltonians including the Born- 
Oppenheimer type Hamiltonian (38), cutting off high energies is equivalent 
to cutting high momenta. Then conservation of energy immediately ensures 
that no momenta exceeding the cutoff are produced over time. This idea 
was developped in [So] and also used in [MaSo]. We will briefly indicate an 
alternative way on how to implement such an energy cutoff in order to fit 
the Born-Oppenheimer and similar settings into our general assumptions. 

Let Hq E S"™ be elliptic and positive, i.e. there is a constant C > such 
that Ho{q,p) > C {p)"^. For example the Born-Oppenheimer Hamiltonian 
as defined in (38) satisfies Hq G Sq and it is elliptic provided that V is 
positive (otherwise just add a constant to Hq since V S S^). Then we can 
prove adiabatic decoupling uniformly for energies below any A G M, i.e. on 
Ran]l(_oo,A](^^o)- 
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Let A = {{q,p) : HQ{q,p) < A}, then bounding the total energy by A 
essentially corresponds to confining the slow degrees of freedom to the region 
A in phase space. More precisely, let xx ^ Co°(^) ^^ch that Xa|[o,A] = 1 
Xa|[a+<5,oo) = for some 6 > 0, then xx{Ho) G 0P5(^°° is a semiclassical 
operator. Furthermore, its symbol x '■= Symb(xA(-ffo)) has an asymptotic 
expansion which is identically equal to 1 on A, i.e. x\a ^ 1 and identically 
equal to on the set where HQ{q,p) > X + 6. The statements about Xxi^o) 
and its symbol follow from the functional calculus for semiclassical operators 
as developped e.g. in [DiSj], Theorem 8.7. 

Next we assume that one can define an auxiliary Hamiltonian -ffaux(<Z)P) £ 
Sq such that 

(i) H^nx{q,p) = Ho{q,p) for all {q,p) G A + S := {{q,p) : Ho{q,p) < \ + S}, 

(ii) H^^^{q,p) > Ho{q',p') for all {q,p) ^A + S and {q',p') eA + S 

(iii) and Ha;ay^{q,p) satisfies the global gap condition (Gap)o. 

This can be easily achieved e.g. in the Born-Oppenheimer setting by replac- 
ing p^ by an appropriate bounded function. 

It follows from the previous discussion that {Hq — i^aux ) XA(-f^o) = 
0_oo(£°°) and that xa(^o) - XA(^aux) = C»-oo(£°°). Using 

Xa(-^o)1(-oo,a](-^o) = ^{-oo,X]{Ho) , 
one finds, in particular, that 

( ^0 - ^aux ) ]l(-oo,A] (^O) = Oo(£°°) (48) 

in the norm of bounded operators and thus also 

(e-^-^^"''*-e-^^°*)]l(_^,,](#o) (49) 

= -ie-iH^u^tj^ d,^iH^u.s (^Q^^ _ ^-iHos ]l^_^^^j(^o) = Oo(£~|t|) . 



Now the scheme of Sections 2, 3 and 4 can be applied to i?aux and by 
virtue of (48) and (49) all results are valid for Hq up to 0{e°°) if one restricts 
to energies below A. In particular one finds that for (g,p) € A the leading 
order symbols of /laux = U*Hg^^^U are given by the formulas obtained in 
Section 4.2 using the symbol Ho{q,p). 
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5 Semiclassical analysis for effective Hamiltonians 



The results of the previous sections are genuine quantum mechanical: semi- 
classical symbols have been used only as a tool in order to construct (and, 
eventually, to approximate) 11 and U, but no semiclassical limit has been 
performed. Indeed, the adiabatic decoupling of energy bands is a purely 
quantum phenomenon, which is, in general, independent from the semiclas- 
sical limit. 

However, under the assumption that CFr{q,p) = {Ej.{q,p)} consists of 
a single eigenvalue of HQ{q,p) of necessarily constant multiplicity i, the 
principal symbol of is a scalar multiple of the identity, i.e. h(){q.p)TTj- = 
E,-{q,p)l!Cf, and a semiclassical analysis of h can be done in a standard 
way. In particular, the dynamics of quantum observables can be approxi- 
mated by quantities constructed using only the classical flow generated 
by the (classical, scalar) Hamiltonian Ej.{q^p). This results in a generalized 
Egorov's theorem, see Theorem 5.1. We emphasize that for more general 
energy bands (Tr{q,p) one cannot expect a simple semiclassical limit, at least 
not in the usual sense. 



5.1 Semiclassical analysis for matrix- valued symbols 

Egorov's Theorem. For the moment, we identify JCf with and h with 
TTr/iTTr, an i X ^-matrix- Valued formal symbol. At least formally, Egorov's 
theorem is obtained through an expansion of the Heisenberg equations of 
motion for semiclassical observables: Let a{q,p,e) G 5°(£, B(C^)), then the 
quantum mechanical time evolution of a is given by 



and satisfies 



da{t) 



-[h, a{t)]. 



(50) 



dt 
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Expanding both sides of (50) on the level of symbols and using [E^l, an{t)] = 
0, 1 = Iqi, one obtains the following hierarchy of equations: 

dao{t) 



dt 



= {E,l,aQ{t)} + i[hi,ao{t)] (51) 



= {EA,ai{t)] + i[hi,a^{t)]-^{{huao{t)}-{a^{t)M}) 

+ i[/i2,ao(i)] (52) 
= Kl,a2(^)} + ^[^l,a2(^)] + ... . (53) 

Since dan{t)/dt does not depend on higher orders, the equations can be 
solved iteratively. The solution of (51) with initial condition 00(^,^,0) = 
aQ{q,p) is given through 

ao{q,p,t) = D*{q,p,t)ao{^\q,p)) D{q,p,t) , (54) 

where : M^*^ — ^ M^'^ is the solution flow corresponding to the scalar Hamil- 

tonian Er{q,p). More precisely, ^*{qo,Po) = {Q{t),p{t)), where {q{t),p{t)) is 
the solution of the classical equations of motion 

q = VpEy. , p = -S/qEy. 

with initial condition (gojj'o)- D{q,p,t) is the solution of 

d 

- D{q,p, t) = -i hi{¥{q,p)) Diq,p, t) . (55) 

with initial condition D{q,p,0) = 1. One can think of (55) for fixed {q,p) € 
M.'^'^ as an equation for the Schrodinger-like unitary evolution induced by 
the time-dependent Hamiltonian hi{^^{q,p)) on the Hilbert space C^. Since 
hi{q,p) is self-adjoint for all {q,p) G M^*^, the solution D{q,p,t) of (55) is 
unitary for all {q,p, t) G M^d ^ 

To see that (54) is indeed the solution of (51), note that the mappings 
U{t) : Cb(M^'^,-B(C^)) ^ Cb(M^'^,-B(C^)) 
defined through (54) for f G M, i.e. 

{U{t) ao){q,p) = D*{q,p,t) ao{^'{q,p)) D{q,p,t) , (56) 

form a one-parameter group of linear automorphisms on the Banach space 
Cb(M2'i, B{0)), since 

{U{s)Uit)ao){q,p) = 

= D*iq,p, s) D*{^%q,p),t) ao{'^' o <f\q,p)) D{^\q,p), t) D{q,p, s) 
= D*{q,p,t + s) ao{^'+'{q,p)) D{q,p,t + s) 
= {U{t + s)ao){q,p). 
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Here the group structures of and of the solutions of (55) are used. Hence 
U{t) is a group and it suffices to check that (54) solves (51) at time t = 0, 
which is easy to see. 

The physical interpretation becomes simpler when translated to states: 
a "classical" particle which started at time at the phase space point {q,p) 
with spinor ipo G C^, is at time t located at the phase space point ^*{q,p) 
with spinor cpt = D{q,p,t)(fo. Hence (55) implies that 

^ = -zh^{<^\q,p))ipt. (57) 

One can also think ofU{t) as being the action on observables of a "classical" 
flow on phase space M^*^ x SU(^) defined as 

mq,p,U) = i^\q,p),D{q,p,t) U) . 



Turning to the higher order corrections (52), (53) etc., they are of the 
form 

= {£^^1, an{t)} + i[hi,an{t)] + /n(ao(t), . . . , a„_i(i)) 

with an inhomogeneity In{t) depending only on the known functions ao(i), 
. . . ,an-i(t). Thus, assuming an(0) = 0, one finds 



«n(*) = / dsU{t - s)In{s) . 

Jo 



(58) 



In order to solve Equation (52) for the subprincipal symbol one needs to 
know However, if one is interested in semiclassical observables with a 
principal symbol which is a scalar multiple of the identity, e.g. in the position 
ao(0) = ql, the last term in ( 52) vanishes at all times, since, according to 
(54), ao{t) is a scalar multiple of the identity for all times. In Section 6 the 
back reaction of the spin of an electron on its translational motion will be 
discussed on the basis of (52). 

We summarize the preceding discussion on Egorov's theorem. 

Theorem 5.1 (Egorov). Let H satisfy either (IG)m for m < 1 and p = I 

or (CO) with p = 0. Let (Triq^p) = {Er{q,p)} be an eigenvalue of HQ{q,p) 

of finite multiplicity i. 

Then the classical flow <I>* generated by Ej.{q,p) and the solution of (55) with 
initial condition D{q,p,0) = 1 exist globally in time. For oq G Sp{B{C^)), 
ao{t) given by (54) is a solution of (51) and ao(i) € 5p(6(C^)) for all t. 
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For each T < oo there is a constant Ct < oo such that for all t G [— T, T] 

\\a{t) - Weiaomi <eCT, (59) 
where a{t) = e*'^*/^ ao e~*'**/^ . 



Proof. Up to the modifications discussed before, the proof fohows easily 
along the lines of Egorov's theorem for scalar valued observables (cf. [Ro, 
BoRo]): To make the expansion of the Heisenberg equation (50) rigorous, 
note that E,- = VTr/iQ'^r G 5*^(1^) with m < 1 and thus the corresponding 
Hamiltonian vector field is smooth and bounded. It follows by standard ODE 
techniques [Ro] that 5tao(^**) € S\ and hence also dtao{t) € 5°, where ao{t) 
is given by (54). Thus one can interchange quantization and differentiation 
with respect to time and obtains 



ait) - We(ao(i)) = f^ds^^ {f''/' We{ao{t - s)) e"^'^^/^) 

^%W,{aQ{t-s))]-We{^{t-s)\ \e 



■ihs/e 



is a semiclassical 



Now, by construction, | h,ao{t — s) —We 

operator in OPS'Ke) with vanishing principal symbol. Hence the integrand 
is really 0(e) as a bounded operator and (59) follows. □ 



This matrix- valued version of Egorov's theorem has been discussed sev- 
eral times in the literature [Iv, BrNo]. 



Berry connection. With this preparation wc explain the motivation be- 
hind the particular splitting of the terms in (37). It is of geometrical origin 
and related to the Berry connection. Recall that in the Born-Oppenheimer 
setting hia/3{q,p) = —p-Aa^{q) and thus Aa/s{q) acts as a gauge potential of 
a connection on the trivial bundle x C^. Its origin is purely geometrical, 
since it comes from the connection which the trivial connection on the trivial 
bundle x H{ induces on the subbundle defined by iroiq). If one assumes 
that Ran7ro(g) is 1-dimensional, the internal rotations along classical trajec- 
tories are just phase changes, the so called Berry phases, and are due to 
parallel transport with respect to the Berry connection [Be, ShWi, Si]. 

In the general case the second term of hi (q, p) in (37) , which we denote 

by 

hBeaf3{q,P) = -i{lpa{q, P), {Er,^0}(.q,p)) , 
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corresponds exactly to this parallel transport along the generalized Berry 
connection. More precisely, the trivial connection on the trivial bundle M^*^ x 
Hi induces a C/(^) -connection on the subbundle defined by 7ro(g,p). After 
unitary rotation uo{q,p) the coefficients of this connection on the bundle 
M2<i X are 

in the sense that a section s{q,p) is parallel if (V — iA)s = 0. It is parallel 
along some curve c(r) = (g(T),p(T)) in M?'^ if 

(dr - c(t) • iA{q{T),p{T))y{q{T),p{T)) =0. 

For classical trajectories, where c{t) = {VpE,., —VqE-^)'^, this condition be- 
comes 

(at + i/iBe(gW,p(t)))s(g(i),p(t)) =0. (60) 

If hi = /iBc, (60) is exactly Equation (57) for the rotation of the spinor 
(pt[q{t),p{t)) = D(^q.,p, t)(p() along the trajectory of the particle. This means 
if hi = hBe, the spin dynamics corresponds to parallel transport with respect 
to the Berry connection along classical trajectories. 

Emmrich and Weinstein [Em We] give a geometric meaning also to the 
remaining terms in their analog of hi . While this is a natural venture in the 
context of geometric WKB approximation, it seems to be less natural in our 
approach, since we work in a fixed basis in order to obtain simple analytic 
expressions. 

Wigner function approach. The previous results on the time-evolution 
of semiclassical observables translate, by the duality expressed through 

('0, aoV') = / Trc<! {ao{q,p)W'^{q,p)) dqdp, 

to the time-evolution of the Wigner transform 

W^iq,p) := Symb(Pv,)(g,p) = (27r)-'^ [ d^e'^'P i^{q + ^riq - 
as 

= / Tr^J ao{q,p)D*{q,p,-t)W'^{^-\q,p))D{q,p,-t)) dqdp + 0{e). 
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Transport equations for matrix-valued Wigncr measures were derived in 
[GMMP] and applied to the Dirac equation in [Sp]. 

Semiclassical propagator. Often one is not only interested in the semi- 
classical propagation of observables, but more directly in a semiclassical 
expansion of the kernel K{x,y,t) of the unitary group 

(e-^'^*/^V)(^)= / dyK%x,y,t)^P{y). (61) 

As in the case of Egorov's theorem, generalizing the known results for 
Hamiltonians with scalar symbols to the case of operator-valued symbols 
is straightforward, whenever the principal symbol of /i is a scalar multi- 
ple of the identity. As in the scalar case, see [Ro], one makes an ansatz of 
the form 

K^{x,y,t) = / dpeH5(-'f.*)-H (^£^a,(x,p,0) , 

where S{x,p,t) is real valued and the aj's take values in the bounded linear 
operators on C^. Demanding (61) at time t = 0, i.e. K'^{x,y,0) = d{x — y), 
imposes the following initial conditions on S and {aj}j>o' 

S{x,p,0) = x ■ p , ao(x,p, 0) = 1 and aj{x,p,0)=0 for j > 1 . 

For later times the coefficients are determined by formally expanding the 
Schrodinger equation for K^{x,y,t) 

ie-^^K'{;y,t) = hK%;y,t) 

in orders of e. At leading order only = Eq contributes and one obtains 
as in the scalar case 

dtS{x,p,t) + E,{x,V:,S{x,p,t)) =0, (62) 

the Hamilton- Jacobi equation for the symbol ho. The next to leading order 
equation is the so called transport equation for oq: 

idtao{x,p,t) = jC{x,p,t)ao{x,p,t) + hi{x,VxS{x,p,t)) ao{x,p,t) . (63) 

The differential operator C{x,p,t) is the same as in the scalar case, sec [Ro] 
for an explicit formula. Here we just want to point out that the known 
techniques from the scalar case apply with one modification: as in (51), also 
in (63) hi contributes as an additional rotation in the transport equation 
for the leading order term. Since the solution of (62) exists only until a 
caustic is reached, the approximation (62), (63) to the propagator is a short 
time result only. The extension to arbitrary times is a complicated task, in 
general [MaFe]. 
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5.2 An Egorov theorem 

Ultimately the goal is to approximate expectation values of observablcs in 
in the original Hilbert space H = L'^{R'^, Hi) rather than in 7i = L^iR"^, Ki). 
Before stating a theorem an obvious, but important observation should be 
made, which seems to have been overlooked, or at least not stressed suffi- 
ciently, in related discussions, e.g., [LiFl, LiWe, BoKe2, MaSo]: We proved 
that in the case CFr{q,p) = {Er{q,p)} the effective Hamiltonian h projected on 
the subspace K, = Ranllr has a semiclassical limit in the sense of a general- 
ized Egorov theorem, in principle, to any order in e. However, the variables q 
and p in the rotated representation are not the canonical variables of the slow 
degrees of freedom in the original problem. More precisely, let qu = x® 
and = —ieVx ® ^Ht be the position and momentum operators of the slow 
degrees of freedom acting on Ti, and let q]Q = and p]Q = —ieVx ® '^K.{ 

be the same operators acting on /C. Then q^c = li^U q^ U* Ilr + 0{e) and 
piC = IlrUpT-cU*Ilr + 0{e), with a, in general, nonvanishing e-correction. 
Physically this means that the quantities which behave like position and 
momentum in the semiclassical limit are only close to the position and mo- 
mentum of the slow degrees of freedom, but not equal. This phenomenon 
is well known in the case of the nonrelativistic limit of the Dirac equation. 
The Newton- Wigner position operator and not the standard position oper- 
ator goes over to the position operator in the Pauli equation. The standard 
position operator has neither a nice nonrelativistic limit nor, as we will see, 
a nice semiclassical limit, because of the Zitterbewegung. Switching to the 
Newton- Wigner position operator corresponds to averaging over the Zitter- 
bewegung, or, in our language, to use the position operator qjc in the rotated 
representation. We remark that in the Born-Oppenheimer case, and more 
generally whenever ttq depends on q only, one has qic = UrU qy, U* nr+0(e^). 

With this warning we exploit that semiclassical observables do not change 
after unitary rotation in leading order and state the Egorov theorem for the 
observables in the original representation. 

Corollary 5.2. Let H satisfy either (IG)m with m < 1 and p = 1 or (CO) 

with p = and let crj.{q,p) = {Ej-{q,p)} consist of a single eigenvalue of 
Ho{q,p) of finite multiplicity i. Let bo £ S^{B{7i^)) such that [bo,TTo] = 
and B{t) := e*^*/*^ 6o e~*^*/^. Let qq := ny^uoboUQ-rrr and define ao{t) is in 
(54)- Then for each T < oo there is a constant Ct < oo such that for all 



For bo = f l-^f , with f G ^^^(M), one obtains as a special case of (64) that 



t e [-T, T] 




(64) 



{B{t) - 6o(^*))n < eCr. 
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Corollary 5.2 follows from Theorem 5.1 and a straightforward expansion 
in e of the terms to be estimated after rotation with U. 



6 The Dirac equation 

6.1 Adiabatic decoupling of electrons and positrons 

We apply the adiabatic perturbation theory to the one-particle Dirac equa- 
tion with slowly varying external potentials, i.e. to 

Hy> = ca- {^-ihVy — -A{£y)^ + Pmc? + e^{ey) 

acting on L^(M^, C^). Here A : — >■ is the vector potential of an external 
magnetic field B = V /\A and : — >■ M the potential of an external electric 
field E = —V(f). For the Dirac matrices a, (3 we make the standard choice 

where a = (ci, (72, (T3) denotes the vector of the Pauli spin matrices. The 
small parameter e > controls the variation of the external potentials. To 
keep track of the size of the error terms, in this section all physical constants, 
including h, are displayed. 

Transforming to the macroscopic space-scale x = ey one obtains the 
Dirac Hamiltonian 

Hd = ca- (^-iehVx - -A{x)^ + [3mc^ + e(j){x) (65) 

and we are interested in the solution of the time-dependent Dirac equation 
for times of order e~^, i.e. in solutions of 

d 

ieh—J)t = HDil)t (66) 
at 

for \t\ = 0{1). The solutions of (66) for small e approximately describe the 
dynamics of electrons, resp. positrons, in weak fields, as in storage rings, 
accelerators, or cloud chambers, for example. 

i?D is the Weyl quantization of the matrix- valued function 
H-D{q,p) = ca- (p- -A{q) \ + /Smc^ + e(f){q) 



G. Panati, H. Spohn, S. Teufel 



185 



on phase space M^, where now Weyl quantization is in the sense of p i-H' 
—ieh^x, i-e. on the right hand side of (78) e must be replaced by eh. h 
appears here for dimensional reasons and is a fixed physical constant. The 
small parameter of the space-adiabatic expansion is e. H-a{q,p) has two 
two-fold degenerate eigenvalues 

E±{q,p) = ±cpo{q,p) + e^(g) 

with the corresponding eigenprojections 



P±iq,p) 



1 ± 



1 



a 



P 



A{q) ) + Pmc 



Po{q,p) 

where po{q,p) = A/rr?c^~+~(p^^"|A(g))2. Obviously 

E+{q,p) - E_{q,p) = 2cpo{q,p) > C{p) > 0, 

whenever A is uniformly bounded. Therefore the corresponding subspaces 
are adiabatically decoupled and the effective dynamics on each of them can 
be computed using our general scheme. Assuming A G C^(M^,M^) and (/> € 
C^(IR^,M), one finds that Hq £ Sl and thus the assumptions from Section 2 
are satisfied. In particular, Hn is essentially self-adjoint on <S(M*^,C^) and 
E± on 5(M^). 

To be consistent with the notation from the previous sections, let TTo{q,p) 
= P+(q,p) be the projector on the electron band. The reference subspace 
for the electrons is /C = L^(M^,C^) and it is convenient to define it as the 
range of 





n, := 



in L^{I 



The only choice left is the one of UQ{q,p) or, equivalently, of a basis 
{V^aiq,p)}a=i,2 of Ran7ro((?,p). Since the degeneracy of Ran7ro(g,p) is related 
to the spin of the electron, a natural choice is the cTx-representation with 
respect to the "mean"-spin S{q,p) which commutes with HY)(q,p) [FoWo, 
Th]. The eigenvectors '(p±{q,p) of the operator 63 • S{q,p) in Ran7ro(g,p) are 



/ (PQ+mc) \ 



Mq,p) = c^f^^, 



+mc) 



V 





V3 

VI +iV2 / 



, ^-{q,p) = c^f^^. 



( \ 

(pQ + me) 



+mc) 



V 



V\ —IV2 

—vz 



J 



We abbreviated p) := c {p— ^A{q)) /pQ [q, p) for the velocity. The relevant 
part of uq for the analysis of the electron band is thus given by u^iq^p) = 
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{il'+{q,p),'tp-{q,p), *, *) with uq G Si- Of course the positron part indicated 
by *'s would be given through charge conjugation. In our construction we 
want to emphasize, however, that no specification is needed in order to 
determine the expansion of the effective electron Hamiltonian he :=llrhllr 
up to arbitrary order. 

An alternative way to arrive at the same uo{q,p) is to note that the Foldy- 
Wouthuysen transformation tiFw(p)) c.f. [FoWo], diagonalizes the free Dirac 
Hamiltonian Hq{p), i.e. ifo with A,(p = 0. Including the fields uo{q,p) = 
upwip— ^A{q)) then diagonalizes H-Q{q,p). 

For the principal symbol of h^ one finds of course 

he,o{q,p) = E+{q,p)lc2 . 

For the subprincipal symbol after a lengthy but straightforward calculation 
our basic formula (37) yields 

Its I Pq 

^Po{q,P) V c{po{q,p) +mc) 

=: -^a-n{q,p). (67) 

Note that the factor h comes from the fact that the rath term in the space- 
adiabatic expansion carries a prefactor /I". Defining 



l{q,p) = 1/a/1 - {v{q,p)/c)'^ =Po{q,p)/{mc) 
one concludes that 

^{q,p) = — B{q) - \ v{q,p) A E{q)) . (68) 

mc \j{q,p) c(l+7(g,p)) J 

We remark that the second term in (37), the "Berry term", does not coin- 
cide with any of the terms in (68). Indeed, the compact expression (68) is 
obtained only through cancellations in more complicated expressions coming 
from both terms contributing in (37). 

We summarize our results on the adiabatic decoupling and the effective 
dynamics for the Dirac equation in the following 



Theorem 6.1. Let A € C^{R-^,R-^) and 4> G C^(M^M). Then there exist 
orthogonal projectors Il± with 11+ + n_ = 1 such that [Hji,Il±\ = Oq{s^), 
and there exists a unitary U and h G OPSj^ with 
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such that 

^-iHnt _ jj* ^-fht u = 0(i{e°°\t\) . (70) 
Here he and hp are semiclassical operators on L^(M^,C^) with 

oo 

he{q,p,s) X E+{q,p)lc2 + ^ £^ K,j{<l,p) 

and he,j = tt,{u# #u*)jTT, G ^^"^'(^(C^)) for all j > 0, where u e S^^{e) 
is constructed as in Section 3. In particular, he,i{q-,p) is given by (67) and 
thus 

he= (c^m^c^ + {-iehV - + e(/)(g)) - £ ^ a • + Oo(£^) ■ 

Analogous results hold for hp. The errors in (70) and (71) are in the norm 
of hounded operators on L^(R^,C^), resp. on L^(M^,C^). 

According to the effective Hamiltonian (71) the (^i-factor of the electron 
equals 2. There would be no problem to add to the Dirac Hamiltonian the 
standard subprincipal symbol [Th], which accounts for the slightly larger 
^-factor of real electrons. Blount [BI2] computes the second order effective 
Hamiltonian /ie,2, which he finds to be proportional to lc2. /ic,2 is a sum of 
terms allowed by dimensional reasoning, i.e. proportional to VS, Vi?, i?^, 
E"^ , EB. Second order corrections seem to be of interest for the dynamics 
of electrons in storage rings. Ignoring the contribution [BI2], nonrigorous 
expansions are [DeKo] and [HeBa]. 



6.2 The semiclassical limit of the Dirac equation 



Equipped with h^fi and h^^i we can apply the general results of Section 5 
on the semiclassical limit to the Dirac equation. Let <&j_ be the Hamiltonian 
flows generated by E±{q,p) on phase space and let B = bl, b e Si(R), be 
a semiclassical observable in the unrotated Hilbert space which does not de- 
pend on spin. From Corollary 5.2 we conclude for each T < 00 the existence 
of a constant Ct such that for all t € [— T, T] 



(s(t)-w,(6($V)) i)nH 



B{t)-We{b{^t))l)u 



< sCt, 



where B{t) = e*^DV{£^) q ^-i-Hiyt/ieh) ^ Hence, to leading order, states in the 
range of H_|_ behave like classical relativistic electrons and states in the range 
of n_ like classical relativistic positrons. We emphasize that, in general, n± 
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are not spectral projections of Hu, since the variation of </> can be larger than 
the mass gap 2m(P. Hence in the limit of slowly varying potentials a natural 
characterization of "electronic" and "positronic" subspaces is obtained which 
does not come from spectral projections of the free or full Dirac Hamiltonian. 

Next we discuss the leading order spin dynamics, which in the first place 
requires to figure out which operator represents the spin of the electron. 
There has been a considerable discussion on this point, cf. [Th], with no 
general consensus reached. We suspect that the problem is void. The wave 
function is spinor valued and what is observed is the spatial splitting of 
different spinor components in inhomogeneous magnetic fields. Hence we 
should pick the "spin observable" S such that the splitting can nicely be 
attributed to it. E.g., in a magnetic field with gradient along the ^-direction 
the eigenvectors of should have the property that their spatial support 
goes either parallel to +z or to —z, but should not split. In view of (71) a 
natural choice is to take as spin operator the vector of Pauli-matrices a in 
the rotated electronic subspace. In the original Hilbert space this amounts 
to 



^ (7 y h 

where S{q,p) is the "mean" spin defined before. 

The leading order semiclassical approximation for 

follows from Theorem 5.1. For each T < oo there is a constant Ct < oc 
such that for t G [-T, T] 

\\a{t)-^)\\< sCt, (71) 
where crojt(g,p, i), k G {1,2,3}, is obtained as the solution of 
d(Tok{q,P,t) i 



dt 2 



a-n{^X{q,p)),aokiq,p,t) (72) 



with initial condition (Tofc((/,p, 0) = (Jk- This follows from the Equations (54) 
and (55) by setting (Jok{q,P,t) = D*{q,p,t) ak D{q,p,t). 

To solve Equation (72) one makes an ansatz crok{Q,P,t) = Sk{q,p,t) ■ a 
with Sk{q,p,0) = efc. Using [(J„,cJm] = '^i^nmk'^k^ one finds that the spin- 
or "magnetization" -vector Sk{q,p,t) is given as the solution of 

^^^^M = -Sk{q,p,t) A n{^\{q,p)). (73) 
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(73) is the BMT-equation [BMT, Ja] on the level of observables. It was 
derived by Bargmann, Michel and Telegdi in 1959 on purely classical grounds 
as the simplest Lorentz invariant equation for the spin dynamics of a classical 
relativistic particle. 

The semiclassical limit of the Dirac equation has been discussed repeat- 
edly and we mention only some recent work. Yajima [Ya] considers time- 
dependent external fields and proves directly a semiclassical expansion for 
the corresponding propagator. As mentioned already at the end of Section 5, 
this program is mathematically rather involved, since one faces the problem 
of caustics in the classical flow, and different expansions have to be glued to- 
gether in order to obtain results valid for all macroscopic times. Based on the 
same approach Bolte and Keppeler [BoKe2] derive a Gutzwillcr type trace 
formula. Since Hd and U* Hd U are isospectral and since (70) holds, a trace 
formula for the eigenvalue statistics of Hj^ could as well be derived from the 
semiclassical propagator of /i = /ic 1 + 1 ^ /ip. As argued in Section 5, the 
latter is somewhat easier to obtain. In [GMMP, Sp] the semiclassical limit 
of the Dirac equation is discussed using matrix-valued Wigner functions. 
Their results hold for an arbitrary macroscopic time interval, but fuse, as 
does the WKB approach, adiabatic and semiclassical limit. No higher order 
corrections seem to be accessible and the results are weaker than ours in the 
sense that the approximations do not hold uniformly in the states. 

This leads us to the next natural question: What can be said about 
higher order corrections? While in general one would need hc^2, according 
to (52) the semiclassical limit of observables of the type b = ^qIc^j ^ 
S'j'(M), can be determined without this explicit information. For such a scalar 
symbol the principal symbol bQ{t), i.e. the solution to (51), will remain scalar 
and thus its commutator with he^2 in (52) vanishes identically for all times. 
The solution bi{t) of (52) with initial condition 6i(0) = 0, is not scalar, in 
general. Hence, at this order there is back reaction of the spin dynamics on 
the translational motion. We illustrate this point for the position operator 
xijliP) = XQ{q,p) := g and refer to [Tea] for a general analysis of the 
higher order effects in the semiclassical dynamics of Dirac particles. Now 
xo{q,p,t) = xo[^\.{q,p)) and xi{t) is obtained, according to Equation (52), 
as the solution of 



with initial condition a:i(0) = 0. The homogeneous part of this equation is 
just the classical translational and spin motion and the inhomogeneity is 



dxi{t) 
dt 



{Ej^l,xi{t)] + i[he^i,xi{t)\ - {/ie,i,xo(t)} 



(74) 



{K,i,XQ{t)} = --a ■ {n,xoit)} , 



(75) 
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which is not scalar and thus responsible for the splitting of trajectories of 
electrons with distinct spin orientation. Hence, as in (58), 



where U{t) is the "classical flow" defined through (56). 

Without claim of rigor, we observe in (68) that for small velocities v(q,p) 
one has 



and thus according to (74), (75) the correction to the velocity is proportional 
to t, corresponding to a constant force with absolute value he/{2mc)\VB\, 



7 Conclusions 

The basic formulae (33), (34) can be applied, in essence in a mechanical 
fashion, to any concrete quantum problem with two provisos. First of all 
the problem has to be cast into the general form (1) and secondly one must 
have sufficient information on the principal sjnnbol Ho{q,p). Depending on 
Hq considerable simplifications of (33), (34) may be in force, one example 
being the effective Hamiltonian of the time-adiabatic theorem studied in 
Section 4.4. As a net result, if the conditions of the space-adiabatic Theorem 
2.1 are satisfied, the full Schrodinger equation is approximated by an effective 
Schrodinger equation referring to a specific relevant energy band. The errors 
are estimated and, in general, the time scale of validity is much larger than 
the one which can be reached within a semiclassical approximation. 

We focused our interest on a single relevant energy band. No information 
on the complement is needed except for global quantities like the resolvent 
{Ho{q,p) — Er{q,p))^^{l — TTo{q,p)). In previous investigations [Bli, LiFl] all 
energy bands are treated simultaneously. An example which would not fall 
under such a scheme is nonrelativistic QED, which governs electrons coupled 
to the quantized radiation field. In this case the principal symbol has a two- 
fold degenerate eigenvalue at the bottom of the spectrum separated by a gap 




n{q,p)^ —B{q). 

mc 

Let us further assume that B{q) = bq^ez-, then 




as expected for a spin-^ particle. 
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from the continuous spectrum, provided \p\ is sufficiently small and there is 
a suitable infrared cutoff [PSTi]. 

The main restriction of our work is the gap condition of Section 2. There 
are two standard mechanisms of how this condition is violated. (1) There are 
two (or possibly more) locally isolated energy bands of constant multiplicity 
which cross on a lower dimensional submanifold. Away from the crossing 
region the wave function in one band is governed by the effective Hamiltonian 
discussed before. If the wave function comes close to the crossing manifold, 
there is a certain probability to make a transition to the other band. In rather 
specific model systems such transitions have been studied in considerable 
detail [Ha2, HaJoi, FcGe, FeLa]. (2) Hq has a smooth band of constant 
multiplicity bordering the continuous spectrum without gap. This is the rule 
in models from nonrelativistic QED with massless photons. Results for the 
massless Nelson model [Tei] indicate that smoothness of iTo{q,p) suffices also 
in general for adiabatic decoupling at leading order with intraband dynamics 
generated by ho + ehi as defined by (36). However, the expansion stops at 
this stage. Physically, the electron looses energy through radiation, which 
means that the next order correction must be dissipative. 

^From the physics point of view the dynamics of molecules and the dy- 
namics of electrons in a solid are the two most prominent areas of application 
for the space-adiabatic perturbation theory. The former has been discussed 
already in Section 4.3. Bloch electrons do not quite fall into our scheme, 
since the classical phase space is R'^ x T'', T'^ a fiat d-dimensional torus. 
This requires substantial changes which are discussed in [PST2]. 



A Operator- valued Weyl calculus 



Pseudodifferential operators with operator-valued symbols have been widely 
discussed in the literature. The results presented in this Appendix can be 
found in [Ho, Fo, Iv, GMS]. We start with some notation. Let £ he a Banach 
space, then C{R^,£) denote the space of £^-valued continuous functions on 
M.'^. In the same spirit we will employ the notation 5(1^"^, £), LP(M.'^, £), with 
the obvious meaning. Note that, in the special case where £ = Tif is an 
Hilbert space, one has ^^(M'^, Hi) = L'^{R'^)®nf. The space of the bounded 
operators on £ will be denoted as B{£). 



192 



Space- Adiabatic Perturbation Theory 



A.l Weyl quantization 

Let A he a ;S(?Yf)-valued rapidly decreasing smooth function on M^'^, i.e. 
A e 5(M^'^, B{TCi)). If we denote by J^A the Fourier transform of A then, by 
Fourier inversion formula, 

where the integral is a Bochner integral for i3('Hf)-valucd functions. This 
suggest to define an operator A G B{7i), called the Weyl quantization of 
A, by substituting e'^^'^^^'^^ with 6*^'''^+^'^) 1-Hf where q is multiplication 
by X and p = —ieVx in L'^{W^). The exponential is defined by using the 
spectral theorem and it is explicitly given by 

(^e'^n-q+^P)^^ (x) = e^<^-^)/2e^^-^V(ic + eO for V e l2(R<^). (76) 



Thus 

A = 



-i-^ / (^A)(r?, (e^(^-^+«-^) l7i,) dr?d^ , (77) 

(27r)" JK2d V J 

and, in particular, 

II^IIb(h) ^ (2^X^, ii(-^^)(^>0IIb(w,) 

which implies that A belongs to B{H) provided the Fourier transform of A 
belongs to B{ni)). We will also use the notation We{A) = A in 

order to emphasize the e -dependence. 

Substituting (76) in (77) one obtains that for every ip G S{W^,'Hi) 

{Ai^){x) = [ A{l{x + y),^) e'^<-yy^i:{y) d^dy, (78) 

i.e. A is an integral operator with kernel 

Taking (78) as a definition, the Weyl quantization can be extended to 
much larger classes of symbols A{q,p). 

Definition A.l. A function A G C°°{M?'^,B{Hf)) belongs to the symbol 



class S'^{B{7ii)) {with m G M and < p < 1) if for every a, /3 G N"* there 
exists a positive constant Ca^p such that 



{d'^,d^,A){q,p) < C^,p (p)-^l' 



sup 

for every p G W^, where (p) = (1 + |pp)^''^. 
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The space S'^iBiJ-Lf)) is a Frcchet space, whose topology can be defined by 
the (directed) family of semi-norms 



f ) = sup sup (p)— {d^d^^A){q,p) 
|a|+|/3|<feq,peMd 



ken. 

B{Hi) 



(79) 

The following result is proved exactly as in the scalar case, cf. also [GMS]. 

Proposition A. 2. Let A G Sp''{B{Hi)), then A given trough (78) maps 
S{R'^,H{) continuously into itself. 



Since A e S'^{B{Hi)) implies A* e ^^(^^(Hf )), the previous result allows 
to extend ^ to a continuous map on >S'(M'^,Hf). 



It is convenient to introduce a special notation for such classes of opera- 
tors acting on S{W^,Hf), called pseudodifferential operators, 

:= {WeiA) : A G 5™(B(7^f))} . 

In the following wc will sometimes denote S^{B{'H{)) simply as and we 
will use the shorthand S"" -.= 3^. Notice that S"^ C 5^ for any p>p'. 

If A belongs to S'^{B{Ti.{)) then the corresponding Weyl quantization is 
a bounded operator on 7f = L^(M'^,7if). The following proposition sharpens 
this statement (see [Fo], Theorem 2.73). 



Notation. Denote by C^{R'^,£) the space of ^'-valued, k times continuously 
differentiable functions on W^, such that all the derivatives up to the order 
k are bounded. Equipped with the norm 



Pile,. := sup sup mA)ix)\\, 



it is a Banach space. 



Proposition A. 3. (Calderon-Vaillancourt) There exists a constant Cd < oo 
such that for every A G C^'^+^{M.'^'^,B{ni)) one has 

ll^lle(^) < Q sup sup {dgd^A){q,p) = Cd\\A\\^2d+i . 



^^^^ |a|+|/3|<2(i+l q,peRd 



B{H{) ""b 



This implies, in particular, that the Weyl quantization, regarded as a map 
We : S^{B{'Hi)) — > B{H), is continuous with respect to the Frechet topology 
on 50(^(Wf)). 
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A.2 The Weyl-Moyal product 

Next we consider the composition of symbols. The behavior of the symbol 
classes with respect to the pointwise product is very simple, as can be proved 
by using the Leibniz rule. 

Proposition A.4. If A e S'^^iBi'Hi)) and B e S^^^iBiHf)), then AB 
belongs to 

S^^+^^{B{Hi)) for every mi,m2 G M. 

The behavior under pointwise inversion is described in the following propo- 
sition. For every T G B{H{) let the internal spectral radius be pint{T) '■= 
inf{|A| : A G a{T)} . 

Proposition A. 5. Assume that A G S^{B{H{)) is a normal symbol which 
is elliptic, in the sense that there exists a constant Cq such that 

Pint{A{q, p)) > Co {pr for any p G M'^ . 

Then the pointwise inverse A~^ exists and belongs to S~"^{B{Hf)). 

Proof. As a consequence of the spectral theorem (for bounded normal 
operators) one has 

\\A~\q,p)\\^^^^^=pUA{q,p))-' <C{p)-"' . 

Similar bounds on derivatives can be obtained by noticing that 

\\^p(^~')\\Bm = \\-A-\V,A) A-Xin,^ < C'ip)—" 
and applying the chain rule. □ 

The crucial result for pseudodiferential calculus is the following. One 
can define an associative product in the space of classical symbols which 
corresponds to the composition of the operators. Given A G {B{7ii)) 
and B G 5'™2(fi(Hf)) we know that A and B map S{M.'^,nf) into itself. 
Then AB is still an operator on (S(M'^, "Hf^ and one can show that there 
exists a unique e-dependent symbol Symb(^B) =:A#Be S^'^+'^^iBiHi)) 
such that 

We{A)Ws{B)=We{A4^ B). 

The symbol A^ B is called the Weyl product (or the twisted product) 
of the symbols A and B. For the proof of the following proposition in the 
operator valued case we refer again to [GMS]. 
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Proposition A.6. Let A G S^^iBiUi)) and B G S'^^{B{ni)), then AB = 
C with C G S'^^+'^^{B{ni)) given through 

C{q,p) = e^p('-^{Vp-V,-V^-Vg)] {A{q,p)Bix,^)) 

=: {A#B){q,p). (80) 

In particular, S')](i3(Wf)) and S^{B{ni)) := (JmGR •^^('^(^f)) algebras 
with respect to the Weyl product #. 

Since the product A # B depends on e by construction, one can expand 
(80) in orders of £. To this end, it is convenient to define suitable classes 
of £-dependent symbols, called semiclassical symbols, which - roughly 
speaking - are close to a power series in e of classical symbols with nicer and 
nicer behavior at infinity. Our definition is a special case of the standard 
ones (see [DiSj, Ma, Fo, Ho]). 

Definition A. 7. A map A : [0,£o) ^ -S'^j^ is called a semiclassical 

symbol of order m and weight p if there exists a sequence {Aj^j^^ with 

Aj G S^~-'^ such that for every n G N one has that (^A^ — X]"=o ^''^j^ 

belongs to uniformly in e, in the following sense: for any A; G N there 

exists a constant Cn,k such that for any e G [0,£o) one has 

n-l 



j=0 

l(m) 



(m—np) 

<Cn,fc£", (81) 

k 



where ||. . .||^ is the k-th Frechet semi-norm in S*™, introduced in (79). 

One calls Aq and Ai the principal symbol and the subprincipal symbol 
of A. The space of semiclassical symbols of order m and weight p will be 
denoted as S'^^e). If condition (81) is fulfilled, one writes 

A^Y^'^A^ in5-(£) 

i>o 

and one says that A is asymptotically equivalent to the series X^j>oS"'^i 
in S^{£). If A is asymptotically equivalent to the series in which Aj = 
for every j G N, wc write A = 0(£°°). To be precise, we should write 
A = 0{s°°) in ^"(e), but the latter specification is omitted whenever it is 
unambiguous from the context. 

In general a formal power series '^j^Q^-' Aj is not convergent, but it is 
always the asymptotic expansion of a (non unique) semiclassical symbol (e.g. 
[Ma]). 
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Proposition A. 8. Let be {AjjjgM cliT' arbitrary sequence such that Aj £ 
S^~-"^. Then there exists A € S'^{e) such that A x X]j>o^"'^i ^^{^) '^'^^ 
A is unique up to 0{£°°), in the sense that the difference of two such symbols 
is 0{e°°) in S'^{e). The semiclassical symbol A is called a resummation 
of the formal symbol X^j>o ^''^j ■ 



The Wcyl product of two semiclassical symbols is again a semiclassical sym- 
bol with an explicit asymptotic expansion (see [Fo], Theorem 2.49). 

Proposition A. 9. If A ^ Ylj>o^''^j '^p^^i^) ^'^^ ^ ^ Yl,j>Q^^^j 
(e), then A ^ B E S'^^~^^'^{e) has an asymptotic expansion given by 

[a # B)^{q,p) = (2i)-'= Yl {i9^9^A,){d^d^Bi)) iq,p) 

\a\ + \f3\+j+l=k ' '"'^'^ 

(82) 

where it is understood that k,j,l G N and a,/? G N'^. 



For example {A # B)o is simply given by the pointwise product AqBo and 

{A # 5)i = ^0^1 + ^1^0 - ^{^0, -Bo} 
where {•, ■} denotes the Poisson bracket on S'^{B{Hi)), defined through 
, , ^ dA dB dA dB 

Notice that, in general, {A,B} ^ —{B,A} since operator-valued derivatives 

do not commute, in particular {A, A} ^ 0. The usual Poisson algebra is 
recovered in the special case in which one of the two arguments is a multiple 
of the identity, i.e. A{z) = a{z)l-}if. 

As a consequence of the previous result, it is convenient to introduce the 
space of the formal power series with coefficients in S^{B{7if)). This space, 
equipped with the associative product given by (82) and with the involution 
defined by taking the adjoint of every coefficient, will be called the algebra 
of formal symbols over B{Hi). In particular we will denote as M™(e) the 
subspace of the formal power series with a resummation in S^{e), i.e. 



j>0 
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In the context of formal power scries, the product defined by (82) will be 
called the Moyal product and denoted simply as 7^ . Notice that # defines 
a map from M^"i(e) x M'^^{e) to M™i+™2(e). The Moyal product can also 
be regarded as a map from M^i(e,;B(Hf)) x Mp{e,ni) to M^^+"'^ {Hi) , 
where in (82) the operator A and its derivatives act on the vector B and its 
derivatives. 

To sum up the previous discussion, we wish to point out that one can 
prove statements on three levels: formal symbols (i.e. formal power series), 
semiclassical symbols, and operators on (S(M'^,7Yf) C L'^(W^,Hi). A sim- 
ple example illustrates the interplay between these levels. Suppose that 
two formal symbols A € MJ^^{e) and B € MJ^'^{e) Moyal commute, i.e. 
[A, B]# = A#B - B#A = 0. Let A^ G S^^{e) and B^ € S^'^{e) be any two 
resummations of A and, respectively, B. Since we know a priori (by Prop. 
A.9) that the Weyl product A^ # B^ belongs to 6'™i+'"2(e) it follows that 
the Weyl commutator [A^, asymptotically close to zero in S^^^"^^ (e), 

which can be rephrased in the following way: for any n, G N there exists a 
constant Cn,k such that for any e £ [0,eo) one has 



[A„B, 



# 



(mi+m2—np) 
k 



If p > we obtain that definitely rrii + m2 —np<0 for some n G N and then 
Prop. A. 3 assures that the operator commutator [74£,i?£] can be bounded 
in the ;B(H)-norm. Moreover, for p > 0, we can conclude that [Ag,^^] 
is a smoothing operator (i.e. it belongs to OPS~°° := HmeROPSJ^) and 
in particular one can prove that it is a "small" bounded operator between 
the Sobolev spaces H'^ and H'^^'^' for any q,r G N. To be precise, for any 
q,r,n eN there exist a constant Cn,q,r such that 



[Ae,B, 



< C £ 



n 



for any e G [0,£o)) where H'^ stands for H'i {W^ ^Tii) . Notice that for p = 
and m\ + m2 =: m > it is not possible to conclude from [A^B]^ = 
that [Ae, Bs\ is a bounded operator, since it could happen - for example - 

that [Air,Bi;\^ = e~^p'^, which is asymptotically close to zero in 5^(e). In 
the following wc will use the same symbol for an element in S'^{e) and its 
expansion in M™(e). As suggested by the preceding discussion, we introduce 
the following synthetic notation. 



Notation. Let be A and B semiclassical symbols in S^{e). We will say 
that B = A + 0-oo{e°^) if 5 — ^4 is asymptotically close to zero in ^"(e) 
for p > 0. 
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With a little abuse, we will employ the same notation for pseudodifferential 
operators too, i.e. we write B = A + C>_oo(£°°) ii B = A + O-ooie"^)^ As 
noticed above this is a strong concept of closeness, since it implies that B — A 
is a smoothing operator. Compare with the following weaker concept. 

Notation. Let be R and S two (e-dependent) operators on V.. We will say 
that R = S-\- Oq{£°°) if for every n G N there exists a constant C„ such that 

11-^ ~ ^\\b(H) ^ ^n^" 

for every e G [0, eo)- In such a case we will say that R is Oq{£°°)-c\osq to S. 
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